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NOTE 

t  This  1.st  Volume  refers  exclusively 

to  the  Primary  Aim  (Theoretical  Study  of 
the  Trajectories;  Mathematical  Model, 

FORTRAN  Program  and  Sets  of  Results)  included 

t  in  the  Second  Year  Proposal  submitted  to  the 

European  Office  of  Aerospace  Research  and 
Development,  in  February,  1983,  that  pro¬ 
duced  G-rant  NO.  APOSR-83-0340  as  a  conti- 

*  nuation  of  the  Grant  NO.  AFOSR-82-0316. 

The  work  done  referring  to  the 
Secondary  Aim  (Icing  Simulator  and  Metering 
System)  is  to  be  presented  in  a  separated 

•  Second  Volume. 


A-4/IMi.9  iMprMiU  tM  INTA 


INTA 


N." 


Pig.  3 


Grant  NO.  AFOSR-8 3-0340 


STUDY  OF  AEROSPACE  MATERIALS,  COAT’INffS,  ADHESIONS 

AND  PROCESSES 


Aircraft  Icing  Processes 
Primary  Aim 


Ernesto  Montiel  R. 

INTA,  Torre  j(5n  de  Ardoz. 
(Madrid) .  Spain. 


14  Seprt ember  1984 


Pinal  Scientific  Report,  15  Sept, 1983  -  14  Sept, 1984 
Airproved  for  public  release;  distribution  unlimited. 


Prepared  for  INSTITUTO  NACIONAL  DE  TECNICA  AEROESPACIAL 
"Esteban  Terradas".  Torrejdn  de  Ardoz,  Madrid,  Spain, 
and 

EUROPEAN  OFFICE  OP  AEROSPACE  RESEARCH  AND  DEVELOPMENT 
London  England. 


▼INI  w  t  ixyi-v 


I  N  T  A 


N. 


P&g.  4 


AGNOWLECrrMENTS 

I  must  express  my  thankfulness  to 
Mr. Viet or  Loranca  Ruiz  and  Mr. Pedro  Lacruz 
Alberquilla  for  their  invaluable  teachings 
and  help  in  using  the  computer  facilities 
of  the  In  Plight  Test  Section,  as  well  as 
to  Mr. Pedro  (Jarefa  Samitier,  chief  of  that 
Section  for  his  authorization  to  use  them. 

I  must  express  also  my  gratitude  to 
my  wife,  for  her  lovely  patience  all  along 
these  two  years  of  work. 


TJ.KI  W  tTM/l-T 


I  N  T  A 


N. 


P^g.  5 


Introduction 

The  main  purpose  conducting  the  work  implied  in 
this  report  is  to  provide  IN'^’A  with  a  set  of  computer 
programs,  included  in  a  coherent  system  for  analyzing 
aircraft  icing  phenomena,  as  a  necessary  help  to  the 
aeronavesability  certification  responsabilities  of 
this  national  Spanish  Institute, 

In  particiilar,  this  report  concerns  to  a  second 
phase  in  which  a  computer  program,  that  provides  the 
trajectories  of  the  droplets  of  a  cloud  when  a  two- 
dimensional  body  of  any  cross  section  moves  normally 
to  its  axis,  is  established  considering  the  case  in 
which  the  Stokes  Law  for  the  motion  of  an  sphere  is 
applicable,  as  well  as  the  case  in  which  that  Taw  is 
not  applicable  due  to  the  higer  He:'molds  numbers 
involved.  In  both  cases  the  assumptions  made  have 
been  the  following: 

-  Low  concentration  of  particles,  in  such  a  way 
thai;  the  flow-field  of  the  air  is  not  affected 
by  the  presence  of  the  particles. 

-  The  effect  of  gravity  is  negligible, 

-  The  effect  of  compressibility  is  negligible, 

-  The  effect  of  the  boundary  layer  is  negligible, 

-  Incompressible,  potential  flowfield.  No  circu¬ 
lation. 

The  trajectories  described  by  the  dropT ets  with 
respect  to  the  section  are  presented,  graphically,  in 
a  normalized  pattern  in  which  the  maximum  dimension 
of  the  section  in  the  drawing  is  always  approximately 
the  same  whatever  the  dimensions  of  the  real  body  may 
be.  Nevertheless,  the  point  where  the  trajectories 
start  can  be  selected  outside  the  enclosing  frame, 
which  is  also  always  the  same. 
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Mathematical  Model 

'T'o  establish  a  M.M.,  corresponding  to  the  trajec¬ 
tories  that  the  droplets  of  a  cloud  describe  with 
respect  to  a  two-dimensional  body  which  approaches 
moving  normalJ-y  to  its  axis,  is  the  first  step-  of 
the  Primary  Aim  established  in  the  November  19811 
INTA  Proposal. 

We  have  established  one  such  Mathematical  Model, 
analyzing  the  relations  existing  between  the  principal 
variables  implied  in  the  physical  phenomenon  of  the 
motion  of  the  droplets  xmder  the  actions  of  the 
viscous  and  inertial  forces. 

The  M.M.  has  been  established  considering  two 
different  cases,  according  to  the  possibility  of 
applying  the  Stokes  Law.  This  possibility  depends 
upon  the  value  of  the  3e3molds  number  that  needs  to 
be  calculated  in  each  one  of  the  steps  in  which  each 
one  of  the  trajectories  is  divided  for  its  calculation 

The  cal CTilat ions  start  from  a  point  that  we  can 
select  arbitrarily  at  such  a  distance  of  the  profile 
that  permits  to  suppose  reasonably  that  the  pertur¬ 
bation  of  the  air  at  this  point,  due  to  the  motion  of 
the  body,  is  neglectable. 

Thou^,  in  both  cases,  we  establish  the  set  of 
equations  mth  a  notation  usual  in  Physics,  afterwards 
in  both  cases,  we  make  a  change  to  an  intermediate 
notation,  by  capital  letters,  as  an  approximation 
to  the  final  notation  to  be  employed  in  the  FORTRAN 
Program  that  constitutes  the  second  step  of  the 
Primary  Aim  of  the  Proposal. 
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Or-.sa  in  •.vAic’i  6Ko\-:es  Ls"i  is 

iili'.n  the  profile  :r.ovo3  in  the  air,  it  oroluces 
in  each  point  of  ohe  fluid  an  induced  v/ind,  which  can 
be  represented  by  a  vector.  The  difference  between 
this  vector  and  the  vector  representing  the  absolute 
velocity  of  the  droplet  considered,  is  another  vector 
U  which  represents  the  relative  velocity  between  the 
air  and  the  droplet.  This  relative  velocity  produces 
on  the  spherical  droplet  of  diameter  "d"  a  force  f 
given  by  the  Stokes  La.w: 

f=3T7^Ud  ( 1 ) 

where  yw  is  the  viscosity  of  the  air.  This  force  has 
the  same  direction  and  sense  that  the  relative  velo¬ 
city.  Then  v/e  can  write: 

f^3yUd  (2) 

This  force,  by  the  other  hand  and  according  to 
the  iTewton's  Second  Law,  produces  an  acceleration  of 
the  droplet  given  by  the  expression: 


f =ra  ■  a 


where  m  is  the  mass  of  the  droplet ; 
m=Volume  •  density=  "3  ^  ['f']  ^ 


Prom  (2), (3)  and  (4)  results: 

3yud=.g-.^j.d3.  r 


U=  a 


when  the  viscosity  y»i  ,  the  density  of  the  droplet 
Pg,  and  the  diameter  “d"  of  the  droplet  can  be 
considered  constants,  the  expression  preceeding  tf 
can  be  represented  by  a  constant  3,  an  then, 
results  the  very  simple  vectorial  expression: 

3-ir  =  r  (6) 

equivalent  to  the  set: 

=  ^x 
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'u  iiiterval  OjT  ':±.r.e  dt 


v/e  :ienor".ir.rte ; 


'"’■•■cf  ’  ^'/-P  “  components  of  the  absolute  velocit  / 

of  the  fluid  at  the  end  of  che  interval. 

Au,.^  ,  ^’^/f  “  at  the  beginnias  of  te  interval. 

u  ^  ,  u  =  comoonents  of  the  a.bsolute  velocity  u 
xg  ’  75  -  g 

of  the  droplet  at  the  end  of  the  interval. 

Au  ,  Au  =  Idem  a‘.  the  beginning  of  the  interval, 

^>3  J  O 

the  component  "x"  of  the  mean  relative  velocity  is 
given  by  the  expression: 

equivalent  to: 

“x  =(*“xf-  2  Au^g+  -  4'S;g)/2  (10) 

because 

U-.,  »  =  AU„  „  +  du  (11) 

xg  xg  xg  ^  ' 

Reasoning  in  a. similar  manner  to  the  "y"  axis 

results : 


“y  *{^V-  "  -'V"  V  - 


with 

"yg  = V  ““yg  <'3) 

The  increment  du._  relates  with  the  acceleration 
xg 

a^  and  the  increment  dt  in  time  by  the  expression: 

^x  ^ 

and  similary: 

a-y  =  du^^dt  (15) 

The  space  traveled  by  the  droplet  along  the  x 
axis  diiring  the  time  dt  cannot  be  expressed  by  the 
simple  formula  dx  =  Au  * dt  because  there  are  zones 

o 

at  which  the  velocity  has  a  very  low  value.  It  is  ne¬ 
cessary  to  use  the  more  complicated  expression: 

=  Au  •  dt  t  ^  •  a  •  dt^  (16) 
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?.n:'.  31"-'  •.-ry: 

.  dt- 

(17) 

Zquc- 1 ions  (  1  o )  " r..i  (17), 

usin;  (lO 

and  (15)  con  be 

presented  os: 

it 

(13) 

dt 

(19) 

or: 
dt  = 

dt  = 


2-dx 


2.AU  -h  du 


2-dy„ 


(20) 

(21) 


2  ^  ^  V 

Substituting  this  expressions  of  dt  in  (14)  and 
(15)  results: 

U,_l/f2.dyJ  (23) 


a 

=  du 

f2-Au 

X 

xg 

V  xg 

a 

=  du,  . 

r  2  •  Au  , 

7 

yg 

V  yg 

By  substitution  of  expressions  (10) 
expression  (7)  we  have: 

B.(au^j  -  2  -""xg+'^xf  -  '^“xg)-  - 

=  d<A^^-(2-Au^^+du^J 

Similary  we  have: 

B-(*V  -  2-Au.^g+Uyj  -  %g)-dyg 

If  we  select  as  direction  of  the  x  axis  that  of 

the  motion  of  the  profile  ,  the  ordenate  of  the 

•nrofile  results  a  constant  and  its  increment  dy^ 

is  always  zero,  then  results  that  the  increment 

dy^  of  the  absolute  ordenate  y^  of  the  droplet 

is  equal  to  the  increment  dy _  of  the  relative 

gn 

ordenate  y^^.^  of  the  droplet  with  respect  to  che 
profile  : 

“  ^gn 


and  (22)  in 

(24) 

(25) 
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In  order  to  obtain  the  air  velocity  components 
due  to  the  motion  of  the  profile  inside  the  cloud, 
in  the  sense  of  the  X  axis,  we  had  proposed  to  use 
the  expressions  given  at  the  next  page  for  the  two 
components  and  .  Those  exoressions  produce 

resul-t^s ,  concerning  non-steady  flow,  that  can  he 
directly  rela+ed  to  the  results  produced  by  the  well 


known  expressions: 

II 

a 

c 

r  1 

1 

CM _ 

1 _ _ 

j  •  COS  9 

(28) 

u^  -  Uo' 

(rf  -  1 

j  •  sin  9 

(29) 

that  provide  the  relative  velocity  components  u 
and  u  in  cylindrical  coordenates  (r,9)  of  the 
air  with  free  stream  velocity  U;-  ,  due  to  the 
presence  of  a  cylinder  of  radius 

The  agreement  between  the  results  produced  by 
these  two  different  methods  is  excellent.  Never¬ 
theless,  when  the  given  expressions  for  and 
are  applied  to  a  sheet  of  zero  thiclcness,  like  the 
case  of  an  infinitely  long  sheet  of  paper,  the  cal¬ 
culated  components  of  the  air  velocity  result  both 
zero,  •'i'his  inconvenience  do  not  appears  when  the 
expressions  of  V  and  are  applied  only  to  the 
segments  that  are  sighted  from  the  point  (x,y)  at 
which  the  velocity  is  being  calculated.  This  is 
the  reason  why  an  "Non  Sighted  Segments  Factor" 

(NS3P)  has  been  introduced  in  the  calculations. 

The  yalue  of  NSSP  can  be  selected  for  each  trajectory 
from  zero  to  one,  and  indicates  de  degree  of  effect 
in  the  results  of  the  segments  that,  though  impli¬ 
cated  in  the  definition  of  the  profile,  are  not 
sighted  from  the  (x,y)  point. 

The  selection  by  the  operator  of  the  value  of 
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NSSP  is  made,  like  the  values  of  other  parameters, 
by  conversation  between  the  operator  and  the  com¬ 
puter. 

The  NSSP  value  is  printed  before  the  initia¬ 
tion  of  each  trajectory,  due  to  its  great  influence 
on  the  results. 

It  is  expected  to  determine  experimentally,  in 
future  research,  the  values  of  NSSP  that  fits  better 
to  the  different  types  of  profiles  (sheets,  cylin¬ 
ders,  airfoils,  etc.)  to  be  considered,  and  perhaps 
the  relative  distance  to  the  profile. 
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rel^-tions  oet'.vesn  gcor'etrical  vari:?.jl;3; 

V  (32) 

=  dy  +  7  (33) 

ss  n  gn. 

As  we  have  said,  if  we  select  as  direction  of 
the  X  accis  that  of  the  motion  of  the  profile  results: 

=  0  (34) 

t  the  space  traveled  along  the  x  aicis 
by  the  profile  during  dt  is  given  by: 

=  (^5) 

As  a  result  of  all  the  preceeding  reasoning,  we 
can  establish  the  following  set  of  equations  as  re¬ 
presentative  of  the  notion  of  the  droplets  of  the 
cloud  produced  by  the  notion  of  the  profile  in  it  : 

3  -(Au^f  -2  Au,^g+  -  du^g  ).  dx  = 


=  du,^g.(2-Au^g  +  du^g  ) 

^■(AUyf  '  2-Au^g+Uy^  -  du^  ).dy 


=-du 

(  2  • 

Au ,  ^ 

■  du 

yg 

yg 

yg 

^xf 

=  '^xf 

’  ^gn 

""yf 

II 

’  ^gn  ' 

Ax 

+  dx„ 

=  dx^ 

+  X 

gn 

sr 

o 

Q  gn 

Ay 

gn 

+  dy 

3 

=  dy, 

^gn 

dx_ 

=  u 

dt 

n 

n 

^^n 

=  0 

dt  = 

2.dx^ 

/(2 

dt  = 

2-d,7^ 

/(2 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


!  is  a 

set  of  10 

equat 

ions  ' 

(ulate 

the  vahacs 

of  the 

10 

^xf 

"V 

'^g 

^n 

du 

yg 

‘iy. 

9 
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.'33ponj.ing  to 

a-n  iiicr  jr.e.it 

the  rela 

tive  ;Ys1 

;anGe  betv*c9n 

tile  droplet  and  tne 

profile 

,  if  the 

values : 

Au 

y  o* 

‘‘“O 

AV 

corresponding  to  the  actual  point  of  the  trajectory 
are  knovm. 

This  set  of  equations  permits  to  calculate,  from 
a  set  of  initial  values,  all  the  trajectory,  step  by 
step,  whereas  the  hypothesis  about  the  validity  of 
the  application  of  the  Stokes  Law  be  correct. 

In  order  to  solve  by  computer  the  proceeding  set 
of  equations,  we  can  introduce  some  transformations 
in  it  ; 

First  of  all  we  shall  introduce  a  new  notation, 
with  capital  letters,  self comprehensive,  by  which  the 
proceeding  set  of  equations  can  be  written  as; 


3*(AUXP-2‘AUXa+U:Q’-  DUXS)OXl  = 


=  OUXC}  V  (  2  •- AU:C1  +  DILIS)  (46 ) 

3‘(AUYF-2*AIjT'1+ irfF-  DUYa)>DY(l  = 

=  DUYd  *  ( 2  « AlTfG  +  DUYG)  ( 47 ) 

UXP  =  UXP  (XGN,YGN)  (48) 

UYF  =  UTP  (XC?r,YaN)  (49) 

AXGN+DXG  =  DXN+XGN  (50) 

AYGN+DYG  =  DYCT+YGN  (51) 

DXN  =  ITfrDT  (52) 

DYN  =  0  (53) 

DT  =  2>DXG  /  (2*-AU:CG  +  DUXG)  (54) 

DT  =  2*DYG  /  (2*AirirG  +  DU'j:G)  (55) 


'  Taking  into  account  that  XC-N  -  AICIIF  =  DXGN  ,  the 

equation  (50)  can  be  presented  as: 

D'll  =  JXG  -  DX3H  (56) 
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-•’roa  (51)  (53)  results  : 

ZH  =  AZ:J:r  +  3Z:-  (57) 

Froui  (52)  results  : 

DT  =  DXN  /UN  (58) 

Proa  (54  )  and  (55 )  results  : 

DYO  =  DX(J*(2*AUYC}+Dirra)/(2*AUXCM-DU;C(J) 

Prom  (54  )  and  (58  )  results  : 

DXN/UN  =  2»DX&/{2''ADXG+DUXCJ)  (60) 

By  substitution  of  DXN  in  the  preceeding 

equation  for  its  expression  obtained  from  (50) 
results  : 

(AXGN  -  XIN)*  (2*AUX{5  +  DUTCG) 


(59) 


DXCJ  = 


2*u?i  -  2»Au:ca  -  Duxa 


Prom  (59)  and  (6l)  results  : 

(AICGIT  -  XGN)*(2-AUYa  +  DUYG) 

jjycj  =  - 

2tUN  -  2*AUXa  -  Dlica 

Ve  will  establish  the  new  variable 
AXGN  -  XC-N 


PACX  = 


2*UW  -  2  AUXG  -  DUZCG 


(61) 


(62) 


(63) 


equivalent  to  ; 


PACX  = 


-DXGN 


2*UN  -  2*AUXG  -  DUZOG 


Prom  (61  )  and  (63-)  results  ; 
DXG  =  ?ACX*(2»AU:CG +DU'XG) 

Prom  (62)  and  (63  )  results  : 

DYG  =  PAG:C»(2'AIJIG  +  DirZG) 


(64) 


(65) 


(66) 
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7:.-o:r.  i;46  ),(47)  (59)  i'ccv.l  .3  : 

-  2'AUT'l^irf7 

DUYT  - -  '  d:j.:9 

allc?  -  2*A'j::'i +  L"i?’ 

We  Gcjn.  establish  the  iiew  variables  : 

3N  =  B»(iLcirr-:ciii)  +  2»un  -  b-auiw 
CP  =  3*(Au:cp-2*Auxc +  uxp)4(A;{;crr-x(iN) 
v/here  : 

AXCN-XCN  =  -DXCN 

Then,  from  (46)  and  (61)  and  taking  into 
account  (68)  and  (69)  results  : 

3N  t  |j] 


(67) 

(68) 

(69) 

(70) 


DIICG  = 


BN^  -  4*CP 


(71) 


We  can  now  solve  the  set  of  equations  by  the 
following  iterative  process,  once  established  the 
value  of  DXCN  : 

-  Calculate  XCN  =  AICIN  +  DXC-N 


-  suppose  a 

value 

of  DYC 

-  Calc\Alate 

YCN  = 

^  AYCN  +  DYC 

-  Calculate 

the  components  UgIP,UYP  of  the 

velocity  at  the 

point  XCN,Y3N. 

-  Calculate 

BN  and  C?  by  (68)  and  (69). 

-  Calculate 

DUXC 

by  (71 ) . 

-  Cal  dilate 

PACX 

by  (64 ) . 

-  Calculate 

D'JYG 

by  (67 ) . 

-  Calculate 

byc 

by  (66 ) , 

-  Compare  the  DY1  c?.lculated  with  the  DY3 
supposed  and  tray  a  nev/  DY3  ,  adequately 
chosen,  till  the  difference  bet'.veen  the  DY'3 
calculated  and  the  JYT  supposed  be  little 
enou_;;h . 

-  Once  detemined  the  correct  value  of  JYG  , 
calculate  the  remaining  variables  : 
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1 


X 


Sat  of  oouaaio.'is  for  bhe 

^enera.1  case 

?oi'  t’.ie  2'3neral  case  (CrC) 

in  v/hicii  Stokes 

La.w 

an  not  necessary  be  applicabl 

e,  we  have  established 

be  following  set  of  equations 

: 

X,  =  Ax  +dx^ 

5r  n‘  cT 

o  a  o 

(72) 

(73) 

•■Si  = 

(74) 

y„  =  Ay  +dy^ 
n  n  n. 

(75) 

=  Ax _ +dx^ 

gn  gn  gn. 

(76) 

^gn  = 

(77) 

^gn  “  ^g  " 

(78) 

^gn  -  ^g  - 

(79) 

dx^  =  dt 

(80) 

dy^  =  0 

(81) 

1  2 

(82) 

iyg  =  Ay^gdt+i.ayg-dt2 

(83) 

(84) 

^“ys  =  ^yg-“ 

(85) 

“xf  =  “xf  '  V  -  V’ 

(86) 

yr  yr  '  gn  ’  •'gn 

(87) 

=((“xf  -  V)  +  [au^- 

1  /  ^ 

(88) 

"y  '  ( ( V  -  V)  +  i^yt  - 

*v) )  /  ^ 

(89) 

^x  =  i  °Dx'7Pf 

(90) 

^y  =  i 

f  =  m  •  a 

X  ‘“g  ‘‘xg 

(92) 

(91) 

f  »  m  •  a. 

y  s  :/s 

(93) 
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=D:-.  -  ''Dx'(-V) 

(94) 

(95) 

R  =  A3S  ( U,  J  d^  /  i/ 

(96) 

Ry  =  ABb  (U,J  d^  / 

(97) 

A  =  n-rf  =  — .  d  ^ 

S  4  s 

(98) 

m 

S  Is  s 

(99) 

The  double  signe  ( 1  )  indicates  that  the  force  acting 
on  the  droplet  is  not  necessary  an  accelerating  force. 
There  are  zones  at  which  the  air  acts  on  the  particle 
like  a  brake  decelerating  its  motion.  It  was  not  ne¬ 
cessary  to  take  this  double  signe  into  account  in  the 
case  of  the  Stokes  Law  because  in  that  case  the  force 
depended  directly  upon  the  relative  velocity  U  due 
to  the  partic^ila^  expression  of  Cjj  versus  Reynolds 
number.  On  the  contrary,  in  the  actual  case,  the  signs 

of  U  and  U  are  not  directly  transmitted  to  the 
X  y 

forces  f  eind  f  because  U  and  U  are  squared 
X  y  X  j 

and  hence  the  necessity  of  the  sign  preceeding  the 
expressions  of  the  forces. 

YVith  respect  to  the  expression  of  the  drag  coef¬ 
ficient,  it  is  different  according  to  the  zone  of  the 
Reynolds  number  at  which  we  are  operating  : 

For  Reynolds  numbers  from  . 1  till  1 .  we  will  use 

the  Oseen's  expression  ; 

24 


+  4.5 


For  Reynolds  numbers  greater  than  1 
our  own  expression  : 


(100) 

.  we  will  use 


10. 


(3.92494-log  Rg)2 
5.322698 


-  .4038228 


( 101 ) 


which  is  valid  for  Reynolds  numbers  from  1.  till  2CC0. 


I  N  T  A 


N. 


Pig.  21 


1 


1 


i’o  3olv2  0/  Gor.puc3r  jhe  331:  of  3quaoic:is  '.ve 
sliall  "lake  3oae  arren^erneats : 

?ir3t  of  all  '.ve  aliall  put  the  set  in  the  inter- 
aelirte  notation  as  follows; 


hi  =  AlCil  -1-  DXd 

(102) 

■rcj  =  AfO  +  DYd 

(103) 

:cr  =  A:or  +  DXN 

(104) 

'O  =  AYir  +  DYN 

(105) 

XltT  =  AXOrr  +  DXGN- 

(106) 

YGN  =  aYGN  +  DYGN 

(107) 

XGN  =  XG  -  XN 

(108) 

YGrr  =  YG  -  YN 

(109) 

DXN  =  UN  *  DT 

(110) 

DYN  =  0. 

(111) 

DXG  =  AUXG  *  DT  +  ACELX  *  DT«»2./2. 

(112) 

DYG  =  AUYG*  DT  -t-ACELY*  Dm2./2. 

(113) 

DUXG  =  ACELX  *  DT 

(114) 

DOYG  =  ACELY  «  LT 

(115) 

UXP  =  DXP  (  XGN,  YGN  ) 

(116) 

UYP  =  UYT  (  XGN,  YGN  ) 

(117) 

UX  =  ((U:CP-UXG)+  (AUXy-AUXG))  /2. 

(118) 

UY  =  ((TYP-UYG)+  (AUYP-AUYG))  /2. 

(119) 

PUERX  =iCDX  *  DEP  *  UX**2.  *  A  /2 

• 

(120) 

PUSHY  =IGDY  »  DEP  ‘  UY*y2.  »  A  /Z 

• 

(121) 

PUSHX  =  LIASbG  »  ACELX 

(122) 

PUSHY  =  LUSbG  k  ACELY 

(123) 

CDX  =  CDX  (RX) 

(124) 

CDY  =  CDY  (RY) 

(125) 

RX  =  ABS  (0X)*DIAMEG  /  VISCI 

(126) 

RY  =  ABS  (UY)»DIA]itEG  /  VISCI 

( 127) 

A  =  PI  k  DIAMEG*k2.  /  4. 

(128) 

MASSG  =  PI  »  DEG  «  DIAMEG*k3.  / 

6. 

(129) 

CD  =  24.  /  R  ^4.5 

(130) 

CD  =  10.fk((3.924g4-AL0G10(R)) 

*k2./8. 

322698 

-  .4038228  (131) 


i 
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in  orlor  t;o  -:ivo 


tlie  fluid  density  in  lb  ft  , 
will  establish,  tlie  coefficient  : 


'.ianoter  in  nicrons  and 


le.sii’ed,  v/e 


C0?U8  =  (DE?I5/l53824.)*(DIAiI(J/lO.)*»  2. 


(132) 


by  which,  the  expression  of  the  force  PX  multiplied 
by  10.  to  8  results; 


PUEEIZ8  =  +  C0Pn8  *  CDX  *  13X»  »  2. 


(133) 


Por  the  Y  accis  7/e  will  obtain  a  similar  expres¬ 
sion  of  PUERY8  : 


PIIERY8  =  +  G0PU8  *  CDY  *  UY*  #  2. 


(134) 


These  variables  PUERX8  .  and  PUEBY8  7/ill  be 
used  to  calculate  the  accelerations  ACELX  and  ..ACELY 
of  the  droplet  along  the  X  and  T  axes  respectively, 
resulting; 

ACELX  =  19.1  ^PUEEIX8  /  (DEG  *  (DIAHG/IO. )  r  *  3.)  (135) 
ACELY  =  19.1  *PUERY8.  /  (DEG  *  (DIAiflG/lO. )  *  *  3.)  (136) 

To  calculate  the  relative  velocity  between  the  air 
and  the  droplet,  we  suppose  that,  for  an  interval  lit- 
tel  enough,  it  is  possible  to  neglect,  in  the  expres¬ 
sions  corresponding  UX  and  UY,  the  differential 
terms  DUXP/2.  and  DDYP/2.  ,  resulting  the  following 
expressions  ; 


UX  =  AUXP  -  AUXG  -  (ADUXG/2.)  •  DTGC  /  ADT 
Ur  =  AUYP  -  AUYG  -  (ADUYG/2,)  »  DTGC  /  ADT 


(137) 

(138) 


Also  would  be  possible  to  use  the  approximate  ex¬ 
pressions  ; 

UX  =  AUX  +  DUX  =  AUX  -*■  ADUX  *  DTGC  /  ADT  ( 139) 

UY  =  AUY  +  DUY  =  AUY  +  ADUY  *  DTGC  /  ADT  (140) 

In  all  these  expressions,  the  fraction  DTGC/ ADT 
takes  care  of  the  first  step  of  this  general  case  (GC) 
for  which  the  value  ADT  of  the  preceeding  DT  has 
not  necessarily  the  value  DTGC  supposed  for  all  the 
steps  in  this  GC. 
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In  the  case  of  FP=1 .  the  rate  of  ice  growing 
on  the  surface  of  the  profile  equals  the  water 
mass  deposition  rate  (VVMDH). 

At  the  zone  determined  by  the  collision  ooints 
(AXCP.AYCP)  *and  (XCP,YCP)  of  two  consecutive  tra¬ 
jectories,  the  Vi/MDR  can  be  calculated,  in  grammes 
per  minute  and  sqiiare  centimeter  by  the  expressions: 

BA33=(  ( AZGP-XCP)*  *  2.  +  {  YCP-AYCP)  2.)*  k  ,5 
mLR=6.*  RLrtC+DVGNOi'  (UN/l 000. )/( BASS  r  100. ) 

where  the  units  to  be  used  are  those  later  indi¬ 
cated  for  each  one  of  the  variables. 

These  t..o  expressions  are  included  as  lines 
0850  and  0851  in  the  AIPO?  Program  to  be  taken 
into  account  in  future  expansions  of  this  program. 
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The  AIP07  FOR'pRAN  _Pro^^ram 

T’his  is  the  program  presented  as  a  result  of 
the  research  work  performed  during  the  period 
15-3ept-83  /  14-aept-84.  It  corresponds  to  the 
second  year  of  work  stated  in  the  1981  Proposal 
concerning  Aircraft  Icing.  All  the  goals  there 
established  have  been  reached. 

The  fundamental  difference  betv^een  the  AIP05 
and  the  AIPO?  programs  lies  in  that  the  former  can 
be  applied  onl.v  to  cylindrical  sections  while  the 
latter  is  applicable  to  sections  of  any  share. 

The  AIP07  program  is  submitted  continuously 
to  a  process  of  revision  in  order  to  introduce  new 
improvements.  This  is  the  reason  why  it  presents 
groups  of  sentences  concerning  future  developments . 

Like  the  AIP05,  the  AIPO?  program  includes  a 
great  deal  of  logical  sentences  constituting  like 
a  matrix  in  which  the  calculating  units  are  immersed, 
giving  to  the  -orogram  a  great  versatility. 

when  3ta~rting  the  program  it  is  possible  to 
select  the  values  of  the  following  variables: 

0=  Chord  of  the  profile  (cm). 

Nb=  Number  of  segments  defining  the  section. 

UNK-  Body  soeed  normal  to  its  infinite  span  (Knots). 
■nSPIS-  Air  density  (lb  ft”-^)  times  10  to  5. 

E'T’AI5=  Viscositv  index  of  the  air  (lb  ft“^  sec”"*) 
times  10  to  5  . 

RL*VC=  LWC  (grammes  meter”^)  . 

DIAMG=  Droplet  diameter  (microns)  . 

DEO-  Droplet  dens it v  (grammes  cm“^)  . 

XONO-  Initial  abscise  of  the  droplet  (cm). 

DV0N0=  For  automatic  case:  Increment  of  the  initial 
ordenate  of  the  droplet  in  order  to  initiate 
a  new  tra.iectory  (cm)  . 
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Q=  For  the  case  in  which  Stokes  Law  is  applicable: 
Number  of  steps  of  value  DZSN  contained  in  the 
initial  abscise  (=  -  XGNO/DX(JN) ,  (adiraensional) . 

COTA  =  Tolerance  in  permitted  to  calculate  the 
solutions  in  the  iterative  processes. 

YGNO  =  Initial  ordenate  of  the  droplet  (cm)  . 

DTGC  =  Differential  of  time  for  the  general  case. 

There  are  other  questions  that,  though  formu¬ 
lated,  are  not  used  in  the  AIP07  Program.  They  are 
put  among  the  initial  questions  to  be  ansy«ered  to 
the  computer  as  a  prevision  for  future  expansions. 
These  questions  are  the  following: 

COMPACT  ICE  DENSITY  IN 

CRAIfflES/CUBIG  CENTIMETER  ?  (DEHIC) 

DENSITY  RATIO  ?  (HDDS) 

iVe  can  answer  1  to  any  of  them  to  pass  to  the 
next  question. 

There  are  other  questions  concerning  to  the  mode 
of  operation: 

It  is  possible  to  calculate  one  trajectory  start¬ 
ing  from  any  arbitrarily  selected  ordenate,  then  an¬ 
other  one,  etc.  ("MANUAL  MODE")  on  condition  that  we 
will  answer  T .  to  the  question  "ONLY  ONE  TRAJEC¬ 
TORY  ?  ",  Also  it  is  possible  to  calculate  automa¬ 
tically  aUthe  trajectories  contained  between  two 
of  initial  ordenates  YG-NO  and  YGNOX  and  seoarated 
DYGNO  at  XGNO  ,  each  one  to  the  next  trajectory,  on 
condition  .that  initially  we  will  answer  0.  to  the 
preceding  question  and  we  will  establish  the  value  of 
DYGNO  when  the  computer  asks  for  this  value.  If 
desired, that  portion  of  the  trajectories  contained 
in  the  reference  frame,  can  be  drawn. 
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The  constant  section  of  the  infinitely  long 
body  to  be  considered  under  the  icing  effect,  can 
be  defined,  at  a  certain  angle  of  attack,  by  the 
abscises  and  ordenates  of  up  to  100  points,  coun¬ 
terclockwise  located  on  the  profile.  This  values 
can  be  given  by  punched  tape,  free  format,  though, 
in  the  case  of  circ^amf erence  or  NACA  four  digits 
profile,  the  program  has  provisions,  for  an  inter¬ 
nal  generation  of  the  defining  points. 

The  AIP07  Program  has  also  provisions  to 
change,  before  starting  the  calculation  of  the 
trajectories,  the  angle  of  attack  in  any  desired 
number  of  degrees  with  respect  to  the  angle  of 
attack  at  which  the  profile  has  been  defined.  If 
desired,  the  ne»»  coordenates  of  the  defining  points 
can  be  given  bv  the  commter  to  the  operator  in 
piinched  tape,  P0RMAT(2X,P1.0.3)  .  3oth,  the  initial 
and  the  final,  angles  of  attack,  as  well  as  the 
change,  are  printed  as  AliP,  ALFA  and  DALP  respec¬ 
tively. 

.The  results  provided  bv  the  computer  consist 

of: 

A) -  The  values  of  the  established  data: 

C , NS , UNK , DEPI5 , ETAI5 , RLWC , DT AMd , DEG , 

DEHIC ,  RDDE ,  XGNO ,  DYGNO ,  Q ,  COTA,  .  YGNO . 

LENS , SLPROM , SLUPTO , COnPLE , DTGC , PROFILE , 

AT.p,ATPA,YGNOX. 

B) -  The  values  of  some  auxiliary  variables,  cal¬ 

culated  in  the  program,  that  can  be  useful  to 
do  some,  always  desirable,  verifications: 

UN  =  speed  of  the  body  (cTn/sec)  . 

VISCI  =  kinematic  viscosity  (stoke)  . 

ETA5  -  viscosity  index  (poise) 
times  10.  to  5. 
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BEP  =  auxij-iarv  variable,  used  for  sizing 
the  drawings. 

SLU  =  auxiliary  variable,  with  no  special 
meaning,  implicated  in  the  Stokes 
case  calculations, 

XCrirTiP  =  left-hand  abscise  of  the  standard 
frame . 

SW  =  signal  of  writing.  It  is  used  for  tabu¬ 
lation  each  SW  steps.  Its  value  is  Q/tOO. 

GOTO  =  per  unit  admisible  error. 

RLAMD  -  auxiliary  variable,  with  no  special 
meaning,  implicated  in  the  calcu¬ 
lations  of  forces  and  accelerations. 

C0PU8  =  Idem. 

DALP  =  Difference  between  the  desired  and 
the  initial  angle  of  attack  of  the 
profile. 

When  we  desire  a  number  of  decimals  greater 
than  the  provided  by  the  condensed  (normal)  table 
of  results,  we  can  obtain  the  values  of  the  same 
variables  under  other  format  with  more  decimals, 
for  the  steps  ranging  between  two  steps,  arbitra¬ 
rily  selected,  of  those  composing  the  tra.iectory, 
on  condition  that  we  answer  T.  to  the  question 
"LENS  ?"  and  give  the  numbers  of  the  initial  and 
the  final  selected  steps  when  the  computer  will 
ask  "PROM  S'^EP  NUMBER  ?"  and  "UP  TO  STEP  NUMBER  ?" 
respectively. 
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G)-  Ghe  values  of  the  varis/oles  at  each  step: 
o  =  a.’,ur.her  of  the  step  in.  the  trajector:/. 

T  =  tiae  in,  seconds. 

XN  =  abscise  of  the  profile  (cm)  . 

X(?  =  abscise  of  the  droplet,  (cm). 

XGN  =  relative  abscise  of  the  droplet 

with  respect  to  the  profile  (cm)  . 

YGN  =  ordenate  of  the  droplet .  ( cm) . 

OX?  =  component  X  of  the  air  velocity. (ciVsec) 
UYP  =  component  T  of  the  air  velocity. ( cn\/sec) 
OXCf  =  component  X  of  the  droplet 's  velocity. 
DYB  =  component  Y"  of  the  droplet 's  velocity. 
DXGr  =  absolute  space  traveled  by  the  droplet 
along  X  axis.  (cm). 

DXBN  =  relative  spabe  traveled  by  the  droplet 
with  respect  to  the  profile  (cm)  . 

DYB  =  absolute  space  traveled  by  the  droplet 
along  r  axis.  (cm). 

REYHM  =  Reynolds  number.  Mean  value  between 
the  end  and  the  beginning  of  the 
actiial  step. 

D) -  The  drawing  of  the  portion  of  the  trajectory, 

or  trajectories,  contained  inside  the  standard 
frame  in  which  the  standard  profile  is 
enclosed. 

E) -  Information  about  \Thether  the  droplet  hits 

or  not  the  profile-  In  the  first  case  the 
legends  "THE  DROPLET  HITS  THE  SURFACE"  is 
printed.  Otherwise  the  legends  "XBNLP 
REACHED"  is  the  one  printed  at  the  end  of 
the  tabxiLation  when  the  droplet  reaches  the 
left  hand  side  of  the  normalized  frame. 
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ra00004  IS  OH  CR00026  USING  00120  BLKS  R«0000 

0001 

FTN4X 

0002 

PHOGRAN  AIP07 

0003 

C 

0004 

C 

AIRCRAFT  ICING  PROCESSES.  RESEARCH  PROCRAH.  INTA. 

0005 

C 

AUrOR:  DR.E.NOHTIEL  R.  ON  BOARD  INSTRUMENTS. 

0006 

C 

PROTECTION  AND  RESCUE  SECTION.  AERODVNAniCS 

0007 

c 

AND  NAVEGABILITY  DEPARTNENT.  IHTA.  SPAIN. 

0008 

c 

THE  A.I.P.  RESEARCH  PROGRAN  IS  SPONSORED  BY 

0008 

c 

THE  AIRFORCE  OFFICE  OF  SCIENTIFIC  RESEARCH. 

0010 

c 

GRANT  NO.  AFOSR-83-0340 . 

001  1 

c 

HITH  THE  COLLABORATION  OF  THE  COMPUTERS 

0012 

c 

LABORATORY  OF  THE  FLIGHT  TEST  SECTION. 

0013 

c 

THIS  PSOGRAFl  CALCULATES  THE  TRAJECTORIES,  OF  THE 

0014 

c 

DROPLETS.  PRODUCED  8Y  THE  MOTION  INSIDE  A  CLOUD 

0013 

c 

OF  AN  INFINITELY  LONG  BODY,  WHICH  CONSTANT  CROSS 

0016 

c 

SECTION  IS  DEFINED,  AT  A  CERTAIN  ANCLE  OF  ATTACK,  BY 

001  7 

c 

THE  ABSCISES  AND  ORDENATES  OF  UP  TO  100  POINTS, 

0013 

c 

COUNTERCLOCKWISE  LOCATED  ON  THE  PROFILE. 

0013 

c 

THIS  POINTS  CAN  BE  GIVEN  BY  PUNCHED  TAPE. 

0020 

c 

THE  LAST  POINT  MUST  BE  COINCIDENT  WITH  THE  FIRST  ONE. 

0  02  1 

c 

FOR  THE  CALCULATIONS  THE  PROFILE  IS  APPROXIMATED  BY 

0022 

c 

A  POLICONAL  LINE  COMPOSED  BY  SEGMENTS  LIMITED  BY 

0023 

c 

THE  GIVEN  POINTS.  BE  AWARE!  THE  NUMBER  OF  THE  POINTS 

0024 

c 

EXCEEDS  ALWAYS  IN  ONE  TO  THAT  OF  THE  SEGMENTS. 

0023 

c 

WHEN  THE  PROFILE  IS  A  CIRCUMFERENCE  OR  A  NACA  FOUR 

0026 

c 

DIGITS,  THE  PROGRAM  HAS  PROVISIONS  FOR  AN  INTERNAL 

0027 

c 

GENERATION  OF  THE  DEFINING  POINTS. 

0028 

c 

IF  DESIRED,  THE  VALUES  OF  THE  MOST  IMPORTANT 

0028 

c 

VARIABLES,  CALCULATED  STEP  BY  STEP  ALONG  EACH 

0030 

c 

TRAJECTORY,  CAN  BE  TABULATED. 

0031 

c 

IF  DESIRED, EACH  TRAJECTORY  CAN  BE  DRAWN. 

0032 

c 

IN  CASE  OF  INPIHGEMENT,  THE  LECEMOE  ’THE  DROPLET 

0033 

c 

HITS  THE  SURFACE.'  IS  PRINTED.  OTHERWISE  THE 

0034 

c 

LEGENDS  AT  THE  END  OF  THE  TABULATION  IS  'XCHLP 

0035 

c 

REACHED. ■ 

0036 

c 

IF  THE  AUTOMATIC  MODE  IS  SELECTED  A  WHOLE 

0037 

c 

SET  OF  TRAJECTORIES  IS  OBTAINED  BETWEEN 

0038 

c 

THE  STATED  LIMITS.  IN  MANUAL  MODE,  EACH 

0038 

c 

TRAJECTORY  MUST  BE  SELECTED  BY  THE  OPERATOR. 

0040 

OINENSION  A(100>,D< 1D0>,XX(190>,YY(I00> 

0041 

C0NN0N/C0H2/X< 1D0>, Y(  1D9 > 

0042 

PI»3  .  141582 

0043 

WRl TE< 12  >0, 1,0.1590 

0044 

WRITE!  12  )1,  1,8888,  1500 

0045 

WRITE(  12  )1,  1,8888, 8500 

0046 

BRITE! 12)1, 1,9,8500 

0047 

WRITE! 12)1, 1,0,1500 

0043 

WRITE! 12 )0, 1 , 0,5000 

0048 

WRITE!l,lO)  - 

0050 

10 

FORHAT!2X,*CYLINOER  ?  ANSWER!  Y£S«1.  M0«0.'> 

0051 

READ(1,«)  ESCYL 

0052 

IF(ESCYL.CT.0.5>  GO  TO  180 

0053 

20 

WRITE! 1, 30) 

0034 

30 

FORMAT!  2X , 'NACA  FOUR  DIGITS  PROFILE  ?', 

0055 

*•  ANSWER  !  YESsl .  H0»0 . '  ) 

0036 

READ! 1 , * )  E$H40 

0057 

IF! ESN4D . GT . 0 . 5)  GO  TO  270 

0058 

40 

NUHBERaO . 

A-4/MttJ  J*!  INTA 


INTA 


N. 
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I 


( 


I 


t 


I 


■  ■ 


0059 

C08T0-0. 

0060 

URirE( 1<50> 

0061 

50 

F0RI18T<2X, -NUMBER  OF  SEGMENTS  ?•> 

0062 

REROCl.*}  HS 

0063 

H-MS*1 

0064 

8RITE< 1,60> 

0065 

60 

F0RM8T< 2S j -rtHGLE  DF  ATTACK  AT  BHICH*/ 

0066 

THE  PROFILE  IS  DEFINED  <DEGREES)  ?■> 

0067 

READ<li*>  ALF 

0068 

HRirE< 1>  70> 

0069 

70 

F0RHAT<2X. -ABSCISES, ORDEMATES  ?  (CM>-> 

0070 

READ<5.* )< A< I >, I-l , 100> 

0071 

READ<5,* )<0<K  >,K-1 , 100) 

0072 

90 

AI-A<  1  ) 

0073 

DO  1 10  J«2,N 

0074 

IF(A<J  )-AI >100, 1 10,  110 

007S 

100 

AI-A<  J  > 

0076 

JAl-J 

0077 

0AI-0<  J  > 

0078 

110 

CONTINUE 

0079 

AO-A<  1  ) 

0080 

DO  130  J*2.H 

0081 

IF<  A<a )-AD>130, 130,  120 

0082 

120 

AO-A<J) 

0083 

JAD«J 

0084 

0AD-0<  J  > 

0085 

130 

CONTINUE 

0086 

0I«0<  1  ) 

0087 

DO  150  J«2,N  - 

0088 

IF<  0< J )-0I >140, 150,  150 

0089 

140 

0I-0<  J  ) 

0090 

JOl-J 

0091 

AOI«A<  J  > 

0092 

ISO 

CONTINUE 

0097 

0S-0<I  ) 

0094 

DO  170  J-2,N 

0095 

IF<0(6  >-03)170,170,  160 

0096 

160 

0S-0< J  > 

0097 

J  OS  *6 

0098 

A0S-A< J) 

0099 

170 

CONTINUE 

0100 

00-<08*0I >72. 

0101 

A0-< AD+AI >72. 

0102 

DIAG-S«RT<< AD-AI >**2 .♦( 03-01 >**2 .  ) 

0103 

C-DIAG 

0104 

DO  180  J>1,N 

0105 

XX( J )-A<  0  >-A0 

0106 

180 

VY<  J )-0<  J >-00 

0107 

GO  TO  470 

0108 

190 

NUNBER-0 . 

0  109 

UR1TE< 1,200) 

0110 

200 

F0RMAT<2X, - INTERNAL  GENERATION  ?- , 

0111 

*■  AN98ER.-  YES-l.  H0-0.-> 

0112 

REA0(1,*>  DIGl 

0117 

IFIDIQl  GT.O  5>  GO  TO  210 

0114 

GO  TO  40 

0115 

210 

URITE( 1,220) 

0116 

220 

F0RHAT<2X, -DIAMETER  (CENTIMETERS)  ?- ) 

0117 

READ<li«>  C 

0118 

CDATO-  1 . 

ttol  IHT4 
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0119 

AOsC/2. 

0120 

3RI rE(  1. 230  > 

0  12  1 

230 

F0R«ftT<2X, ■HOUSER  OF  SECHENTS  ?  ■, 

0122 

•  ■<LEOTER  THAN  100>'  > 

0123 

REAO< 1 .«  >  NS 

0124 

N-HS*'! 

0125 

ZET1»2 

i 

0126 

A( 1 >»A0 

0127 

0<1)«0. 

0123 

J=1 

0129 

00  240  I-l.NS 

0130 

J  z  J  4.  1 

i 

0131 

A< J  >aAO«COS(ZETI*l  > 

1 

0132 

240 

0< J  )sA0*SIN(2ETI*l  > 

0133 

URITE<  1. 250  ) 

0134 

250 

F0RHAT<2X, "ANGLE  OF  ATTACK  AT  NHICH’ 

0135 

«■  THE  CYLINDER  IS  DEFINED  (DECREES) 

?■  > 

0136 

READ<1>*>  ALF 

0137 

DO  260  4=1, N 

0133 

XX<  J  )  =  A( 4  ) 

0139 

260 

YYC  4  )  =  0<  4  ) 

0140 

CO  TO  470 

0141 

270 

yRITEC  1,230) 

, 

0142 

230 

F0RHAT<2X, ■ INTERNAL  GENERATION  ?■  , 

0143 

•  ■  ANS«ER :  Y£S=1  .  H0  =  0 . ■  > 

1 

0  144 

REAO( 1 , «  )  D  IS2 

0143 

1F(0IG2.CT.0.5)  GO  TO  290 

0  146 

CO  TO  40 

0147 

290 

URirE<  1,300  ) 

i 

0143 

300 

F0RHAT<2X . 'NUNBER  OF  THE  PROFILE  ?*, 

0  149 

«■  (FOR  EXANPLE:  2409>’> 

0  130 

READd,*)  H'JBBER 

0151 

310 

WRITEC 1,320) 

0  152 

320 

F0RHAT(2X, 'CHORD  ?  (CH)*) 

0133 

REA0(1.«)  C 

1 

0134 

COArO«l . 

0135 

330 

VRITEC  1,340  > 

0156 

340 

FORRAKZX, 'ANCLE  OF  ATTACK  AT  BHTCH 

THE', 7, 

0157 

*■  PROFILE  IS  DEFINED  (DECREES)  ?  . 

ANSHER:  0.') 

0133 

READd. «)  ALF 

0159 

L1*IFIX(NUNBER/1000.  > 

0160 

K3»NUI18ER-1000*L1 

0161 

L2»IFIX(  K3/IOO.  > 

0162 

K2aK3-100*L2 

0163 

YA=Ll*C/100. 

0164 

XAaLZ^C/IO. 

0165 

THICK»K2*C/100. 

0166 

DELTX»C/20. 

0167 

X«-DeLTX 

0163 

DO  430  4*1,50 

0169 

El»(3./4.-l ./200. >*C-X 

0170 

IF( El . CT . 0)  CO  TO  350 

0171 

E2«(  7. /8 . -1 . 7400. >*C-X 

0172 

IF(E2.GT.O)  CO  TO  360 

0173 

E3*(  13.716. -1  .7800.  )*C-X 

0174 

1F(E3.CT.0)  CO  TO  370 

0175 

E4»(  31 .732. -1  . 71600  .  >*C-X 

0176 

IF(e4.CT .0)  CO  TO  380 

0177 

E5*( 63. 764.-1 .73200. >*C-X 

0178 

1F(E5.CT.0)  CO  TO  390 
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0179 

0  180 

0181 

0182 

0183 

350 

0184 

0185 

360 

0186 

0187 

370 

0188 

0'89 

380 

0190 

0191 

390 

0192 

0193 

400 

0194 

410 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

420 

0  20  2 

0  20  3 

0204 

0205 

0206 

430 

0207 

0208 

0209 

0210 

440 

0211 

0212 

0213 

0214 

0215 

450 

0216 

0217 

0218 

460 

0219 

0220 

0221 

0222 

0223 

461 

0224 

0225 

0226 

0227 

0228 

470 

0229 

480 

0230 

023  1 

0232 

0233 

490 

0234 

0235 

0236 

500 

0237 

0238 

E6»<  127.  <^128.  -  I  .  >*C-X 

IF( E6. GT . 0>  GO  TO  400 
OELTX-C/1  152 
GO  TO  410 
DELTX-C/20. 

GO  TO  410 
DELTX-C740. 

GO  TO  410 
OELTX-C/80. 

GO  TO  410 
OELTX-C/160  . 

GO  TO  410 
DELTX-C/320. 

GO  TO  410 
OELTX-C/640 . 

X-X+DELTX 

BET-C-X 

C0-86S<8ET/C> 

TERCO-1 .4845*SQRT(C0> 

YT«THICK*<  TERCO-CO*< . 63  +  C0*(  1 . 758-C0*( 1 . 421 5-CO* . 5075  >  > >  > 
AYC=YC 

IF<  X8-8ET )420. 420,  430 
YC-Y8*X*<  2. *<  C-X8>-X)/(C-XA>**2 . 

TZ«<  YC-ftYC  )/l>ELTX 
C0SZ«1  ./'SQRT(  I  .+T2»*2.  > 

SIHZ»T2*C0S2 
GO  TO  440 

YC«YA*<C-X)*(  2 .•XA-C*X>/XA**2  . 

T2-< YC-ftYC)/OELTX 
COSZ-l ./SQRT( 1 .+TZ**2.  ) 

SIHZ  — TZ>fCOSZ  ‘ 

CALL  XYUL<X, YC.YT, SIHZ, COSZ, XU, YU , XL , YL) 

A<a  >-XL 
0<4>-YL 
I-lOl-J 
A<I  )-XU 
0<1>-YU 

DO  480  J-1, 100 
XX<  4  )-A<  J  )-C/2  . 

YY(  J  )-0(  J  ) 

ALF-0. 

H-100. 

HS-H-1  . 

HRITE<  1,461) 

F0RnAT<2X,*HUST  I  PRINT  THE  ABSCISES  AH0”,7,- 
OROEHATES  OF  THE  PROFILE  AT  THE  ANGLE*,/, 

OF  ATTACK  AT  UHICH  IT  HAS  BEEN  DEFINED  ?',/) 

READ<1,*>  PAD 
IF<PA0.LT.0.5)  GO  TO  530 
URITE<6,4dO)  ALF 

F0RHAT(lHl,/'/,2X, 'ANGLE  OF  ATTACK  AT  IIHICH', 

*’  THE  PROFILE  IS  DEFINED  )',X,F6.2,X, 

♦■DECREES.',///) 

HRITE<  6,  490  ) 

F0RHAT<4X, 'DEFINITION  OF  THE  PROFILE', 

*■  AT  THIS  ANCLE  OF  ATTACK:',/) 

HRITE<6,S00) 

FORHATI IHO, 4X, 'ABSCISES  OF  THE  POINTS', 

*■  LIHITIHC  THE  SEGHEHTS : ' , / > 

«RITE<  6,  520  >  A 
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0239 

URITE<6.510> 

0240 

510 

FORMATS IHO, 4X, 'OROENATES  OF  THE  POINTS", 

0241 

LIHITIHG  THE  SEGMENTS > 

0242 

UR1TE<  6,  520  )  0 

0243 

520 

F0RMAT(:5<2X.F8.2>> 

0244 

530 

yRlTE<l,  540  > 

0243 

540 

F0RMAT<2X . "OESIREO  ANGLE  OF  ATTACK  ?", 

0  24  4 

*■  (DEGREES)") 

0247 

REA0(1.4)  ALFA 

0249 

DALF-ALFA-ALF  , 

0249 

0ALFR-0ALF*PI/180.  1 

0250 

DO  550  J>1,H 

0251 

X< J  >-XX<  J )*C0S(DALFR)-YY< J  >*SIH<DALFR> 

0252 

550 

Y(J )«XX< J )*SIN(DALFR)+YY( J )*C0S(DALFR> 

0253 

J  *1 

0254 

X-X<  9  ) 

0255 

Y-Y( J> 

0256 

CALL  ORAU<C,X,y,NXP,HYP  ) 

0257 

HR1TE< 12)0. l.NXP.HYP 

0259 

DO  560  J>2.H 

0259 

X-X( J  > 

0260 

Y=Y<  J  ) 

026  1 

CALL  DRAy<C,X. Y.NXP.NYP  ) 

0262 

560 

yRITE<  12  )1. l.NXP.HYP 

0263 

yRITE< 12 )-l . 1.9999.5000 

0264 

570 

X0»X<1  ) 

0265 

DO  590  J-2.M 

0266 

IF<X<J )-XD>590.390.580 

0267 

580 

XD-X<J) 

0269 

JXD«9 

0269 

YXD-Y<J) 

0270 

590 

CONTINUE 

027  1 

DO  591  J-l.H 

0272 

X< J  )«X< J  )-XD 

0273 

591 

Y< J >-Y<J >-YXD 

0274 

yRIT^<1.592) 

0273 

592 

F0RHAT<2X."«0ST  I  PRINT  THE  ABSCISES  AND". 7. 

0276 

ORDEHATES  OF  THE  PROFILE.  AT  THE  DESIRED"./. 

0277 

*"  ANGLE  OF  ATTACK.  HITH  RESPECT  TO  THE"./. 

0279 

*■  ACTUAL  AXES  ?"./) 

0279 

REA0<1,«>  PXY 

0290 

IFCPXY.lt. 0.5)  GO  TO  600 

0281 

HRITE<6.593>  ALFA 

0292 

593 

F0RHAT<IH1,//.2X. "ACTUAL  ANCLE  OF  ATTACKi". 

0283 

*X,F6. 2.x. "DEGREES.",///) 

0294 

HRITE<  6.  594  > 

0285 

594 

F0RNAT<4X. "DEFINITION  OF  THE  PROFILE". 

0286 

*•  AT  THIS  ANGLE  OF  ATTACKi"./) 

0287 

yRITE<6.S9S) 

0289 

595 

FORHATC 1H0,4X, "ABSCISES  OF  THE  POINTS". 

0299 

*•  LIMITING  THE  SEGMENTS:"./) 

0290 

yRITE<  6.  597  >  X 

0291 

URITE<6.596) 

0292 

596 

F0RNATC1H0.4X. "ORDEHATES  OF  THE  POINTS". 

0293 

*"  LIMITING  THE  SEGNEHTS : ". / ) 

0294 

yRITE<  6,  597  )  Y 

0295 

597 

F0RMAT<5(2X.F8.2>> 

0296 

600 

XI-XC 1 > 

0  29  7 

DO  602  J-2.N 

0299 

IFCXC J  )>XI )  601.602,602 
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0299 

601 

XI«X<J) 

0300 

JXI>J 

0301 

0XI>0< J) 

0302 

602 

CONTINUE 

0303 

610 

URITE<  1.620  ) 

0304 

620 

F0R«8T<2X, 'BODY  SPEED  IN  KNOTS  ?* > 

0305 

RE8D<li*>  UNK 

0306 

URITEC  1.630  ) 

0307 

630 

F0RHrtT<2X. 'rtlR  DENSITY  IN  LBCHl/'CUBTC  FEET*. 

0308 

*’  TINES  10  TO  5  ?• > 

0309 

RE80< 1 .•  )  0EF15 

0310 

URITE(  1.640  ) 

031  1 

640 

F0RH«T<2X, 'VISCOSITY  INDEX  IN  LB( H  )/< FT*SECONO  )• . 

0312 

•*  TINES  10  TO  5  ?• > 

0313 

RE8D< 1 .«  )  ET8I5 

03M 

URl TE<  1. 650  ) 

0315 

650 

F0RnAT<2X. "LIQUID  M8TER  CONTENT  IN', 

0316 

•'  GRRHNES/CUBIC  METER  ?•) 

0317 

REA0<1.*)  RLUC 

0313 

URI TEC  1. 660  > 

0319 

66  0 

FORMATC 2X DROPLET  DIAMETER  IN  MICRONS  ?' ) 

0320 

REAOCl.*)  DIAMG 

0321 

BRITE<  1.  670  ) 

0322 

670 

F0RMAT<2X, 'DROPLET  DENSITY  IN', 

0323 

*'  GRAMMES/'CUBIC  CENTIMETER  ?') 

0324 

REAOC I .*  )  DEG 

032< 

URITEC  1.680  ) 

032« 

630 

F0RMAT<2X. 'COMPACT  ICE  DENSITY  IN', 

0  32  7 

•'  GRAMH6S/CUBIC  CENTIMETER  ?') 

0328 

READ< 1 ,♦  )  DEHIC 

0329 

URITE<  1.690  )  ! 

0330 

690 

F0RHAT<2X, 'DENSITY  RATIO  ?') 

0331 

READCl.*)  RODE 

0332 

yRITE<  1.  700  ) 

0333 

700 

F0RHAT<2X, ' INITIAL  ABSCISE  OF  THE  DROPLET  NITH', 

0  33  4 

*/.•  RESPECT  TO  THE  BODY  IN.  CENTIMETERS  ?',/. 

0335 

•'  RECOMMENDED:',/. 

0336 

*'  XGNO-IO*C*) 

0337 

READCl.*)  XGNO 

0338 

URITEC  1.710) 

0339 

710 

F0RHAT<2X, 'FOR  AUTOMATIC:  INCREMENT  IN', 

0340 

♦'  CENTIMETERS'./.'  OF  THE  INITIAL  ORDENATE'. 

0341 

•'  OF  THE  DROPLET'. /.'  IN  ORDER  TO*. 

0342 

•'  INITIATE  A  NEB  TRAJECTORY  ?') 

0343 

READCl.*)  DYGHO 

0344 

URITEC  1. 720  > 

0345 

720 

FORHATC2X,'FOR  STOKES:  ABSC XGNO/DXGN  ) 

0346 

•'  CIQUAL  OR  GREATER  THAN  100.)') 

0347 

READC 1 ,*  )  Q 

0348 

BRITEC 1.730) 

0349 

730 

F0RHATC2X, 'PER  CENT  TOLERANCE  IN  THE  EXPLORATION  ?' ) 

0350 

READCl,*)  COTA 

0351 

yRlTE<  6.  740  ) 

0352 

740 

F0RMATCIH1.2X.  'DATA:'.//) 

0353 

BRITEC  6.  750  ) 

0354 

750 

F0RHATC6X, 'C . 7X,'HS' ,7X, 'UNK', 5K. 'DEF15' , 

0355 

*4X.  'ETAIS',  5X,  'RLBC'.SX,  'DIAMG*. 5X,' DEC, 

0356 

*5X.  'DEHIC,  4X.  'RDDE*.  6K.  'XGNO*.  4X,  'DYGHO*  , 

0357 

•6X, 'Q' ,8X, 'COTA' ,/ > 

1  0359 

IFC CDATO . LT .0 . 5)  GO  TO  770 

I  N  T  A 


N. 


P4fl.  13 


0359  yRITEC  6. 7  60 )C. NS<UNK.DEFI5<ETAI5>  RLVC, OIAHG/ 

0  36  0  tiOEC,  OEHJ  C  ,  RDOE  ,  XGHO  i  OYGNO,  a,  COTA 

0  36  1  76  0  F0RI1AT<4X,F5.2/5X.F4.0.4X.F4.0.5XiF5.0,4X, 

0362  **F5.2.5X.F5.3.5X.F4.1.4X>F4.3,5X,F4.3.5X, 

0363  ♦F4.3.6X,F40.5X.F3.1.5K.F6.0.3X,F7.6.7/) 

0364  GO  TO  790 

0365  770  yRITE<6. 780)HS.UHX,DEFI5.ETAI5/RLUC/DIAHG/ 

0366  ♦OEG- OEHIC-ROOE.XGHO.OYQMO.Q.COTA 

0367  780  FORMATC 1 4X; F4 . 0. 4K, F4 . 0. 5Xi F5 .0i4X. 

0363  *F5.2.5X.F5.3>5X,F4.1,4X,F4.3.5X/F4.3,5X, 

0363  *F4  3 ,6X. F4 . 0, 5X.F3 . I, 5X, F6 .0, 3XiF7 .6,// ) 

0370  790  URITEUiSOO) 

0371  800  F0RHAT<2X. ‘ONLY  OHE  TRAJECTORY  ?  (HAHUAL  HDDE)'. 

0372  *■  AHSyER!  YES*1.  .  H0=»0«> 

0373  REAOd.*)  SUT 

0374  AUTOM-l.-SUT 

0375  CN-0. 

0376  yRITE<1.808> 

0377  808  F0R«AT<2X, 'FOR  AUTOKATICt  HAXIHUH  VALUE*./, 

0378  *•  ALLOyED  FOR  THE  INITIAL  ORDENATE  ?* > 

0379  REAOd,*)  YGHOX 

0330  809  URITEd.dlO) 

0381  810  F0RH^r<2X, 'INITIAL  OROENATE*. 

0332  ••*  OF  THE  OROPLET  IH  CEHTIHETEftS  ?*> 

0383  R£A0<1.*)  YGHO 

0384  IF< YGH0>811.818.312 

0385  3U  T0P£»-<C  +  YX1>> 

0336  IFIYGHO.GT.TOPE)  GO  TO  813 

0387  IF<SUT.GT.0.5)  GO  TO  813 

0388  GO  TO  813 

0339  812  TOPE-C-YXO 

0390  IF( YGHO. GT. TOPE)  GO  TO  815 

0391  CO  TO  818 

0392  81  5  yRlTECliSU) 

0393  816  F0RHAT<2X, 'ATEHTION  !  YGHO  OUT  OF  THE*./, 

0394  ♦*  HORIZONTAL  LIMITS  OF  THE  GRAPHIC*./, 

0395  *•  00  YOU  HAHT  PRINTED  RESULTS  ?*,/, 

0396  *•  ANSWER!  YES*1.  .  H0»0.*./> 

0397  REAOd.*)  PR 

0398  IFCPR.GT .0.5)  GO  TO  818 

0399  GO  TO  809 

0400  818  HRITE(t.819> 

0401  819  F0RHAT<2X,*FAST  CALCULATION  FOR  REYNH  LEATER  THAN  .1  ?*./. 

0402  •*  CAOTIONi*./, 

0403  «*  ANSWER  1.  IS.  IH  GENERAL.  NOT  RECONHENDEO . * . /. 

0404  *•  ANSWER!  YES-1.  .  NO-0.*) 

0405  READ<1,*>  SUD 

0406  URITE<  1.820  ) 

0407  920  F0RI1AT<2X,*LEHS  ?  YES-l.  ,  NO-0.*) 

0408  R£A0<1.*)  LENS 

0409  HRITE<1.830) 

0410  830  F0RMAT<2X, *FRON  STEP  HUMBER  ?*> 

041 1  READ4 1 >«  )  SLFROM 

0412  HRirE<1.340) 

0413  840  F0RMAT<2X, 'UP  TO  STEP  HUMBER  ?*  ) 

0414  READ<1.«>  SLUPTO 

0415  WRITEd.SSO) 

0416  850  F0RMAT<2X, 'TABULATION  IN  COUPLE  ?*, 

0417  •'  YES-l.  .  H0-0.*> 

0418  REA0<1.«>  COUPLE 


lMt»r«uU  lUl  INTA 


I  N  T  A 


N. 


Pia.  36 


0419 

VRITEC  1. 860  > 

0420 

960 

F0R«Ar<2X, ‘DT  FOR  THE  GENERAL  CASE 

0421 

♦'  RECOMHEHDEDi‘,7. 

0422 

*•  OrGC*< XGNO/20000  )/UNK  •/ / > 

0423 

READ< 1 . ♦ >  OTGC 

0424 

IF< NUH8ER . LT. 0 .5 >  GO  TO  890 

0425 

8RI T£<  6/  870  > 

042% 

870 

F0RMAT<5X/‘SUT*/5X,'YGH0'/6X,*LEHS*/5X/ 

0427 

f'SLFftOM*  ,  6X,  'SLUPTO*  /  6X/*C0LIPLE*/7X/’DTGC'/ 

0428 

♦  6X,  'PROF ILE •/ 8X, 'ALF* , 7X/ "ALFA* / 6X, -YGNOX*/ /  ) 

0429 

HR1TE(  6..  880  >SUT,YGN0/LENS/SLFR0f1/SLUPT0, 

0430 

«COUPLE. DTGC/ NUMBER / ALF/ ALFA/ YGHOX 

0431 

880 

F0R«AT<6X/F2. 0/5X/ F5  1/6X/F2 .0/ 6X/F5 .0/ 

0432 

*7X/F5.0/9X/F2.0/8X,F6.5/8X/I4<6X/ 

0433 

•  F5.  1 /6X/F5  1/6X/F5.  1,//) 

0434 

GO  TO  920 

0435 

890 

URITE<6/900> 

0  438 

900 

FQRHAT<5X. ’SUT' , 5X/ 'YGHO * , 6X / ‘LEHS'/SX/ 

0437 

•  •SLFR0ri'/6X/'SLUPT0*/6X/  'COUPLE*/ 7X/  'OTGC*/ 

0438 

*7X/  * ALF" / 3X- *  ALFA* / 6X , * YGHOX * , 7  > 

0439 

HRI TE<  6/ 910 )SUT/ YGHO/ LENS/ SLFROM/ SLUPTO, 

0440 

♦COUPLE /OTGC/ ALF/ ALFA/ YGHOX 

044  1 

910 

FORKAT<  6X , F2 . 0 , 5X/ F5 . 1/ 6X, F2 . 0/ 6X/ F5 .0/ 7X/ 

0442 

♦  F5  .  0 , 9X/  F2 . 0  /  8X/ F6 . 5/ 6X/ F5 . 1 / 6X/ F5 . 1 / 6X/F5. 1 / 7/> 

0443 

920 

WRI TE<  6/  9  30  > 

0444 

93  0 

F08MAT<2X, “AUXILIAR  PARAMETERS  CALCULATED  BY*, 

0445 

*•  AIP07  PROGRAHi •,//> 

0448 

UH«5l , 48*UHK 

0447 

VISCI-929 . *ETAI5/&EFI5 

0448 

ETA5-I4. 895*ETAI5 

0449 

9EP-C 

0  45  0 

DXGH--XGN07Q 

0451 

GLU--1800  0  .*&XGH’*ETA57(  DEG*D  I  AnG**2  .  > 

0452 

XGHLP--<  4 .5*C/7  +  XD  > 

0453 

TOPOS-C-YXD 

0454 

SW-Q/lOO .  ■ 

0455 

COTO-COTA7IOO . 

0458 

RLAHO-OI AHG/( 10000 . *VISC I > 

0457 

C0FU8-< DE FI  571 5832  4 .  >*(  D rAHG7 10  .  >**2 . 

0458 

URIT£<  6/  940  > 

0459 

940 

F0RHAT<6X,*UH*/7X/ "VISCI*/5X/  'ETAS'/GX, 

0460 

•  •BEP*/ 6X/ 'GLU* ,6X, 'XCNL?*, 6X,  ’SU* , 8X. 

0461 

••GOTO* /7X,*RLAH0*/ 6X, ‘COFUB* / 9X,*DALF*, 7) 

0462 

URI TE<  6/  950  )UH. VISCI/ ETA5/ BEP, 

0463 

♦GLU/ XGHLP ,SU/ GOTO/ RLAMD , C0FU8 , DALF 

0464 

950 

F0R«AT<4X,r7.2,2X/F7.4/4X/F6.2/3X/ 

0465 

*F6 . 2/ 3X/ F7. 2/ 3X/ F6 . 2/ 4X, F4 . 0/ 5K/ 

0466 

•F8. 7,4X/ F7. 6/ 4X/F7 . 6. 6K/ F5 . 1  I  77) 

0467 

URITE< 1/960) 

0468 

960 

F0RHAT<2X.*NSSF  ?  <071.  )*> 

0469 

READ<l/«)  CHSS 

0470 

969 

URI TE<  1/  970  ) 

047  1 

970 

F0RHAT<2X, ‘TABULATION  ?  !  YES-1.  ,  NO-0.*) 

0472 

read<1/4')  tab 

0473 

artab-tab 

0474 

980 

«RITE<l/  990  ) 

0475 

990 

F0RNAT<2X. ‘TRAJECTORY  IK  PLOTTER  ?*, 

0476 

•  ■  ANSHER:  YES-1.  /  NO-0 .  * > 

0477 

READ<l/*>  PLOT 

0478 

HRI  TE<  1/  1.000) 

*-t/UiRS.S  ImprMU  «l*l  INTA 


INTA 


N.* 


Pig. 


S7 


0479 
0430 
048  1 
0482 
0483 
0484 
0485 
0486 
0487 
0483 
0489 
0490 
0491 
0492 
0493 
0  49  4 
0495 
0496 
0497 
0498 
0499 
0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 
0508 
0509 
0510 
051  1 
0512 
0513 
0514 
0515 
0516 
0517 
0518 
0519 
0520 
0521 
0522 
0523 
0524 
0525 
0526 
0527 
0528 
0  52  9 
0530 
0531 
0532 
0533 
0534 
0535 
0536 
0537 
0538 


1000  F0RM8T<2X, 'DO  YOU  M8HT  8,0  IM*, 

*•  PUNCHED  TAPE  ?  :  YES»1.  ,  HO-0.*) 

REAi><l,«>  PUNCH 

I F( PUNCH . LT . 0 . 5  )  GO  TO  1 020 

tfRI TE<  4,  t  010  )( A<  I  >,  I-l ,  1  00 ),  <  0( K  >, K-1, 100  ) 

1010  FORaAT<2X,F10  3> 

1020  CONTINUE 

1030  YGNO-YGNO-OYGNO 

1040  S=0. 

AYGNO» YGNO 
YGHO-YGHO+OYGNO 
IF<  YGHO-TOPOS  >1050, 1049,1049 
10  49  IF(  YGN-TOPQS>1050, 2290,  2  290 

1050  IF<YGNO-YGHOS >10  52,  1052,  1051 

1051  I F< AUTON . GT . 0 . 5 >  GO  TO  2310 

1052  U«0. 

DT-OTGC 
T-0  . 

XN-0  . 

OXC-0 . 

XG-XGNO 

OYG-0. 

YG»YCHO 

DXCN--XGNO/Q 

XGH>XGNO 

OYGN-0. 

YGH-YG 

CALL  SIG<:NS,XGN,Y&N,CN3  3,3IG1,SIG2,SIG3,SIG4> 

VX«<  SIGl  'SIG3  )4>UH/PI 
VY«<  SIG2-SIG4 )*UH/PI 
UXr »VX 
UYF* VY 
DUXG>0 . 

UXG-UXF 

OUYG-O. 

UYQ-OYF 

ARX-0. 

ARY-0. 

AREYM-0 . 

RX-0  . 

RY-0  . 

REYN-0 . 

REYHfl-0. 

SF-0  . 

FT«0  . 

TSUP-3. 

ACORTE-CORTE 

AXCP-XCP 

AYCP-YCP 

ABKCP-8XCP 

ABYCP-BYCP 

CORTE-0. 

1060  HRITE<6, 1070>  YGHO,CHSS 

1070  F0RMAT</7,4X, * YCHO-', X, F4 . I, 6X, 'NSSF*' , X, F5 . 3, // > 

URITE(6, 1080) 

1080  F0RMAT<6X,’S’,8X,’T'.7X, ■XN',7X,"XC’,6X, 

*‘XGN*,  6X,  ‘YGM*  .6X,  ■UXF"  ,  6X.  'UYF’,  6X,  ‘UXC  ,6X, 
**UYG‘, 6X, 'OXG* ,6X, *DXGH* ,6X, ■DYG" ,6X,*REYNH'  ,  /  ) 
1090  1F(  TAB . LT . 0 . 5 >  GO  TO  1110 

URI TE< 6.  1  too >  3.T,XN,XG,XGN, YGH,UXF,UYF, 


A-4AMI2.S  lmpr«uU  titti  INTA 


INTA 


N 


Pig.  38 


OSH? 

♦UXG.UyG.OXG.OXGH.OYG, REYHM 

054') 

1 1  00 

F0Rn8T<2X,F6.0,4X.F7.5.2X/F6.1/2X/F7.5.2X, 

054  1 

*F6.1.3X.F7.4,X,F8.2iX,F8.2,X,F8.2,X/ 

0542 

*F8.2,2X,F7.5.2X,F7.2/2X,F8.5/2X,F9.3) 

0543 

1 1  10 

S-S  +  1  . 

0544 

ORX  =  RX-A.OX 

0545 

ORY«RY-ARY 

054i4 

ORE VN-REYN-AREYH 

0547 

AOT=DT 

0548 

AT-T 

0549 

AXH-XH 

0550 

AXG«XG 

0551 

ADYG-OYG 

0552 

AYG*YG 

0553 

AXGH-XGN 

0554 

AYGN-AYG 

0555 

AUXF-UXF 

0556 

AUYF-UYF 

0557 

ADUXG-OUXG 

0558 

AUXG-UXG 

0559 

AOUYG-DUYG 

0560 

AUYG-UYG 

0561 

ARXoRX 

0562 

ARY-RY 

0563 

AREYH-REYH 

0564 

V>H>1 . 

0565 

IF<SM-f.2-W>  1  1  20,1  130,1  1  30 

0566 

1120 

«»1  . 

0567 

IF<S.LT. 1.5)  GO  TO  1630 

0568 

IF< ORX  .  LT . 0 . 0 . OR . ORX . EQ . 0 . 0 >  GO  TO  1630 

0569 

IF( DRY  LT  .0 .0 . OR . ORY. EQ . 0 . 0)  GO  TO  1630 

0570 

I F< OREYM . LT . 0 . 0 . OR  .  DREYN  .  EQ . 0 . 0  )  GO  TO  1630 

0571 

1130 

IF<AREYH-.l)  1139,1139,1630 

0572 

1139 

IF<  SUO .LT  .0 .5  >  CO  TO  1630 

0573 

1140 

OXGH  — XGHOXQ 

0574 

XGH-AXGH+OXGH 

0575 

IFC XGH-XGNLP>  1180,1190,1150 

0576 

1150 

SZC-O. 

0577 

IFC  XGN  )  1  160,  1  160,  1290 

0578 

1160 

IFC  XI-XCN  )  1170,  1170,  1290 

0579 

1170 

SZC-l . 

0580 

GO  TO  1290 

0581 

1180 

XGH-XGMLP 

0502 

1190 

IFC  SF. LT , 0, 5)  GO  TO  1260 

0583 

1200 

URITEC  6, 1210) 

0584 

1210 

FORflATC/,  7X,*XGNLP  REACHED*,//) 

0585 

MRITEC 12>-1, 1,9999, 5000 

0586 

FT-1  . 

0587 

CORTE-0. 

0588 

IFC  AUTOH . CT . 0 . 5>  GO  TO  2250 

0589 

IFC  OPT .GT ,0 . 5  >  GO  TO  1240 

0590 

DPT-l . 

0591 

GO  TO  1240 

0592 

1211 

VRITEC 1, 1212) 

0593 

1212 

F0RHATC2X,*HSSF  ?*,/> 

0594 

READC1,'»>  CNSS 

0595 

1220 

URITEC 1, 1230) 

0596 

1230 

F0RMATC2X. ’YCHO  ?*,/> 

0597 

REAOCl,*)  YGHO 

0599 

IFC  YGH0)123l,  1233,  1  232 

k-kfUULl  laprMiU  INT4 


1  N  T  A 

□[ 

N.*  P4a.  39 

0599 

1231 

T0PE=-<C+YX0> 

0600 

IF( YGNO. GT. TOPE)  GO  TO  1233 

0601 

IF<  SUT .GT .0.5)  GO  TO  2299 

0602 

GO  TO  1233 

0603 

1232 

TOPEaC-YXD 

0604 

IF< YGNO. GT. TOPE)  CO  TO  2299 

0605 

1233 

WRITE! 1, 1235) 

0606 

1235 

F0R(18T<2X,  ■REMEMBER!  YOU  CAN  DE  TERMI  NE  ■  ,  /  > 

0607 

*■  A  TAHGENCY  POINT  BY  LOOKING  EXCLUSIVELY’./, 

0603 

*•  AT  THE  VALUES  GIVEN  BY  THE  PRINTER  FOR’./. 

0609 

*’  THE  COLLISION  POINTS  CORRESPONDING  TO’./. 

0610 

*•  DIFFERENT  VALUES  OF  SOME  YGNO  ADEQUATELY’,/. 

0611 

*•  SELECTED  IN  MANUAL  RODE.  THEN  :’./. 

0612 

*’  IS  THIS  TRAJECTORY  EXCLUSIVELY  TO’./. 

0613 

♦’  DETERMINE  A  TAHGENCY  POINT  ?’./, 

0614 

*’  ANSWER!  YES  =  1  .  .  N0  =  0.’./, 

0615 

*■  (NO  MANDATORY,  ONLY  RECOMMENDED)’./) 

0616 

READ!!.*)  TTOTP 

0617 

IF(TTDTP.LT.0.5)  CO  TO  969 

0618 

PL0T*O . 

0619 

TAB=0. 

0620 

PUHCH=0 . 

0621 

CO  TO  1030 

0622 

1240 

WRITE!  1.  1  250) 

0623 

1250 

F0RMAT!2X, ’ANOTHER  TRAJECTORY  WITH  NEW  MSSF’./. 

0624 

*’  AND/OR  NEW  YGNO  ?  ANSWER:  YES»1.  ,  N0=0.’,/) 

0625 

READ!1.*)  ATWNY 

0626 

IF(ATWMY.GT.0.5)  CO  TO  1211 

0627 

CO  TO  1270 

0623 

1260 

SF*1. 

0629 

CO  TO  1290 

0630 

1270 

WRITE! 1, 1280) 

0  631 

1280 

F0RHAT<2X , ’ANOTHER  TRAJECTORY  WITH  HEW  DATA  ?’ . 

0632 

•’  YES=1.  N0=0.’,/) 

0633 

REA0<1.«)ATHHD 

0634 

IF!  . 5.LT . ATWND  )  CO  TO  530 

0635 

CO  TO  2320  . 

0636 

1290 

DYG=ADYG+AOUYG*AOT 

0637 

CALL  SUDYG! DYG . XGN . AYCN , NS . UN , DXCN , AUXG . AUYG . 

0638 

*AUXF. AUYF.CLU, CNSS. COTA, COYC, DUXC. FACX. DUYC. OXF,  UYF ) 

0639 

OYC=COYG 

0640 

OYCN*OYG 

0641 

0XGsFACX*<2 . * AUXG* DUYC ) 

0642 

XG-AXG*DXG 

0643 

YC=AYC+DYC 

0644 

YCM=YC 

0645 

DXH»OXG-OXCN 

0646 

XN>AXN>DXN 

0647 

DT>DXN/UN 

0643 

T=AT*OT 

0649 

UXGsAUXODUXC 

0650 

UYC*AOYC+DUYC 

0651 

U»! ! UXF-UXC  )**2  . ♦! UYF-UYC  )**2 . >** . 5 

0652 

REYH»U*DI ARC/! lOOOO .•VISCl > 

0653 

REYMH«<REYH*AREYN)/2. 

0654 

1300 

IF! .S-S2C)  1310, 1310, 1301 

0655 

1301 

IF! SF- .5  )1820,  1820,  1811 

0656 

1310 

IF! TSUP. LT. 0. 5)  CO  TO  1390 

0657 

IF<TSUP.GT.2.  )  CO  TO  1370 

0658 

1320 

J>J2*1 

viNi  w  t-smft-v 


I  N  T  A 


Pio-  40 


0659 

1330 

! 

XUEP-X<J) 

0660 

IF(XGH-XUEP>  1340, 1350; 1360 

066  1 

1340 

J  ■  J  +  1 

0662 

OHJ-H-J 

0663 

IF<0HJ .GT .0.5>  GO  TO  1330  ; 

0664 

J*1 

0665 

GO  TO  1330 

0666 

1350 

YIJEP-Y<J> 

0  66  7 

IF(YGN-YUEP)  1520, 1460; 1470 

0668 

1360 

XUJl-XUEP 

0669 

J  1  -  J  ; 

0670 

YUJ1-V<J> 

0671 

j.j-1  1 

0672 

XUJ2-X( J  > 

0673 

YUJ2-Y< J ) 

0674 

J2-J 

0675 

GO  TO  1480 

0676 

1370 

IF( YGH>  1380, 1440, 1450 

0677 

1380 

TSOP-0. 

0673 

J2-JX0 

0679 

1390 

J-J2-1 

0660 

1400 

XLEP-XCJ) 

0681 

1F<XGH-XLEP)  1410,1420,1430 

0682 

1410 

J-J-1 

0683 

IF<  J .CT. 0 . 5  )  GO  TO  1400 

0684 

0685 

GO  TO  1400 

0686 

1420 

YLEP-Y<J) 

0687 

IF(YGH-YLEP>  1470,1460,1540 

0688 

1430 

XL4 1-XLEP 

0689 

J1»J 

0690 

YLJ 1-Y<J  ) 

0691 

J*4  +  l 

0692 

XLJ2«X<J>  1 

0693 

YLJ2-V<J>  ! 

0694 

42-J  ' 

•0695 

GO  TO  1500 

0696 

1440 

XGH-0. 

0697 

GO  TO  1460 

0698 

1450 

TSUP-i. 

0699 

J2>JXD 

0700 

GO  TO  1320  ! 

0701 

1460 

GO  TO  1560 

0702 

1470 

GO  TO  1820 

0703 

1480 

1F(XU42.EQ.XUJ1>  GO  TO  1820 

0704 

TERM5»< YUJ2-YUJ1 >/< XUJ2-XUJ1 > 

0705 

EXGN«rERH5*<XGH-XUJ 1 >+YUJl 

0706 

IF<YGH-EXGH)  1490,1460,1820 

0707 

1490 

COIL  CORTO<  XUJ  1, YUJ 1, XU4  2, YU J 2, AXGN, OYGN, 

0708 

*XGH, YGN,XCP,YCP> 

0709 

GO  TO  1570 

0710 

1500 

IF<XL42.EQ.XLJ1>  GO  TO  1820 

0711 

TER«5-( YLJ2-YLJ1  )7<  XL J2-XL J 1 > 

0712 

EXGH-TER«5*<XCH-XLJ 1 >+YLJl 

0713 

IF(YGH-EXGH>  1920,1460,1510 

0714 

1510 

COLL  CORTO(XLJ  1. YLJ 1, XLJ2, YL42.0XGN,0YGH, 

0713 

*XGH, VGH, XCP, YCP> 

0716 

CO  TO  1570 

0717 

1520 

4H1 -J-l 

0713 

IF( /HI .GT .0 .5  >  60  TO  1530 

A-4/IWS.)  «!•!  INTA 


INTA 


N 
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0713 

1 

J=M 

0720 

1530 

J=J-1 

0721 

XUEP2»X<  J  ) 

0722 

r’JEP2  =  Y<  J  ) 

0723 

COLL  C0RTA( XUEP, YOEP. XUEP2i YUEP2. flXCN. AYCN, 

0724 

•  XGM , YCH, XCP , YCP  > 

0725 

CO  TO  1570 

0726 

1540 

DNJ*M-4 

0727 

1F( OHJ  CTOS)  CO  TO  1550 

0728 

J=1 

0729 

1550 

0730 

XLEP2aX(J) 

0731 

YLEP2aYC J  ) 

0732 

CALL  CORTACXLEP, YLEP,XLEP2. YLEP2. AXCN. AYGN, 

0733 

*XCN, YCN,XCP,YCP> 

0734 

CO  TO  1570 

0735 

1560 

XCP^XCN  i 

0736 

YCP=YCN 

0737 

BXCPsAXCN 

0738 

BYCP=AYCN 

0739 

1570 

XCX=XCP 

0740 

VGH=YCP 

074  1 

IF(AUT0«-.5)  1580.1580^1980 

0742 

1580 

IF<PL0T.LT.0.5)  GO  TO  1610 

0743 

CALL  0HAy2<  C. XO, YXO.XCH. YCH, HXP.HYP) 

0744 

1590 

URITEt 12)1. l.NXP.HYP 

0745 

1600 

yRI  TEC  12 )-l .  1 . 9999, 5000 

0746 

1610 

UR1TE<6. 1620) 

0747 

1620 

F0R«AT<7.6X,*THE  DROPLET  HITS  THE  SURFACE. ■,/> 

0748 

TAS^O. 

0749 

FTal  . 

0750 

C0RTE=1 . 

075  1 

GO  TO  2060 

0752 

1630 

OT*OTGC 

0753 

UX»AOXF-AUXC-C AOUXG/2 . >*OTCC/AOT 

0754 

UY»AUYF-AUYG-<  AI)UYC/2  .  )*DTCC/AOT 

0755 

U*SQRT<UX**2.+UY**2. > 

0756 

REYN=RLAI1D*U 

0757 

REYN««(REYN+AREYN>/2. 

0758 

RX>RLAnD*ASS( UX) 

0759 

RY»RLAHD*ABS(UY) 

0760 

IF<  RX)1658.  1658.  1637 

0761 

1637 

IF( . I-RX>1639. 1639. 1638 

0762 

1638 

C0X=24./RX 

0763 

GO  TO  1660 

0764 

1639 

IF< 1  .-RX  >1650.  1650,  1640 

0765 

1640 

C0X»24 . /RX>4 . 5 

0766 

GO  TO  1660 

0767 

1650 

CDX«10 .•«<(  3. 92494-AL0C10( RX  >  >««2 . /8 . 322698- . 4038228 > 

0763 

1658 

IF<  RY>1690,  1690,  1660 

0769 

1660 

1F(  .  1-RY  >1662.  1662.  1661 

0770 

1661 

CDY«24 ./RY 

0771 

GO  TO  1690 

0772 

1662 

IF(  1  .-RY  >1630, 1630,  1670 

0773 

1670 

C0Y«24.7RY*-4.5 

0774 

GO  TO  1690 

0775 

1680 

COYalO.««<(3. 92494-AL0C10<  RY  >  )**2 . /8 . 322698- . 4038228 > 

0776 

1690 

1F( UX>  1700,  1705.  1710 

0777 

1700 

F0ERX8»-C0FU8*C0X*UX**2 . 

0778 

GO  T.O  1720 

la^ttuU  iM  INTA 


INTA 

N.* 
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d779 

1705 

1 

1 

FUERX8-0.  ■  ■' 

07B0 

GO  TO  1720 

078! 

1710 

FUERX8-C0FU8*C0X*UX**2. 

0782 

1720 

IF(  UY>1730,  1735.  1740 

0783 

1730 

FUERY9«-C0FU8*CDY*UY**2  . 

0784 

GO  TO  1750 

0785 

1735 

FUERY8-0  . 

0786 

GO  TO  1750 

0787 

1740 

FUERY8>C0FU8*CDY*UY**2 . 

0788 

1750 

ACELX-19  .  H>F0ERX8/<0EG*<:DIAHG/'10.  >**3.  ) 

0789 

ftC£LY«19.  l*FUERY3/(OEG*(OIAHG/10.  >**3.  ) 

0790 

0UXG-8CELX*DT 

0791 

DIJYG-8CELY*0T  i 

0792 

0XG><AUXG>DiJXG/2.  >*DT 

0793 

0YG*< AUYG+DUYG/2 . >*0T 

0794 

0XN*UH*0T 

0795 

OXGH-OXG-DXN 

0796 

OYGN-OYG 

0797 

XGH>8XGH>DXGH  i 

0793 

YGH-AYGH+OYGH 

0799 

CALL  SIG< NS. XGH. YGN. CHSS . SIGi .SIG2.SIC3 ,SIG4  ) 

0800 

VX»< Siai-SIG3 >*UH7PI 

0901 

VV-< SIG2-SIG4 >*UH/Pr 

0802 

LtXF=iVX 

0803 

UYF-VY 

0804 

XG-8XG>DXG 

0805 

YG-8YC+0YG 

0806 

XN-8XH+0XN 

0807 

T-RT+OT 

0808 

UXG>8UXG^0UXG 

0809 

OYS-ftUYG-fOOYG 

0810 

1760 

I F(XGN-XGHLP> 1790. 1800. 1770 

081  1 

1770 

IF<  XGN )1780. 1780. 1820 

0812 

1780 

I F(  XI -XGN >1970, 1970, 1820 

0813 

1790 

XGH-XGHLP 

0814 

1900 

IF< SF. LT . 0. 5)  GO  TO  1810 

0815 

GO  TO  1200 

0816 

1810 

SF-1. 

0817 

1811 

XC-XGNLP 

0818 

YC-YGH+( XGHLP-XGN>*DYGH/DX6H 

0819 

XGH-XC 

0820 

YGH-YC 

0821 

1820 

IF(PLOT.LT.0.5)  GO  TO  2040 

0822 

C8LL  0R8«2<  C. XD. YXO ,XGN . YCH. NXP. HYP> 

0823 

CALL  0R8W2( C. XD. YXD. AXGN.PYGH.HXPA.HYPR  ) 

0824 

IF<NXP. LT. 9999. OR. NXP. EQ. 9999)  GO  TO  1830 

0925 

8RI TE( 12 >-I , 1 , 9999. HYP 

0826 

GO  TO  2040 

0827 

1930 

I  F< HYP -0500 >1840. 1850. I860 

0828 

1840 

I  F<  HYP8- 8500  >1870..  1940,  1910 

0829 

1850 

IF<  H  YPrt-8 500  >1940,  1950..  1950 

0830 

1860 

IF(  NYPA-8500>1910,  1950,  1950 

0831 

1870 

1F( N VP-1  5 00 >1890, 1890, 1900 

0832 

1890 

IF<HYPA-1500>1960, 1860,1920 

0833 

1990 

I F(  HYP  A- 1500 >1960, I960,  194  0 

0974 

1900 

IF(  HYP  A- 1500 >1920,  1940,  1940 

0835 

19  10 

COE FY-<HY PA- 8500  )7<  8500 -HYP > 

0836 

HYP-8500 

0837 

GO  TO  1930 

0836 

1920 

C0EFY-<HYPA-l50O>/'<  1500-HYP> 

A-V0«t3  lM|>r«uU  (M  INTA 
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0839 

HVP»1500  1 

0844 

1930 

HXP»<COEFY*HXP  +  NXPA  )/<  DOEFY*-!  .  > 

084  1 

1940 

URl TE< 12  >l . 1 . NXP .HYP 

0842 

GO  TO  2040 

0  84  3 

1950 

URITE<  i2)-l.l. NXP. 8500 

0844 

GO  TO  2040 

0845 

1960 

URI TEC  12  >-1.1. NXP. 1500 

0846 

GO  TO  2040 

0847 

1970 

GO  TO  1310 

0848 

1930 

CN-CN+1. 

0849 

IFC  1  .5-CN)1990. 1990. 1580 

0850 

1990 

B8SE«<  <AXCP-XCP>**2 .♦( YCP-AYCP>**2 . . 5 

0851 

C 

HNDRa6*RLyC«l>YGH0«C  UH/IOOO  .  > 7C  8 ASE*  1 00  .  ) 

0852 

IFC YCP.GT.O.O)  GO  TO  2020 

08S3 

2000 

GO  TO  1530 

0854 

2020 

GO  TO  1580 

0855 

2040 

IFC  2 .2-«  >11  10.  11  10.2050 

0856 

2050 

IFC COUPLER  1 .2-H>ll  10.  1  1  10.  20  6  0 

0857 

2060 

IFC  TAB .LT .0 .5  >  GO  TO  2100 

0858 

URITEC6.2070>S.T.XH.XG.XGN.YCN.UXF,UYF, 

0859 

*UXG , UYG. DXG .DXGH. OYG, REYNM 

0860 

2070 

F0R«ATC2;<,F6 . 0 ,4X.  F7 . 5. 2X.  F6  .  1 . 2X.  F7 . 5. 2X. 

086! 

*F6.1.3X.F7.4.X,F8.2,X.F8.2.X.F8.2.X. 

0862 

*F3. 2.2X. F7. 5. 2X. F7 . 2. 2X. F8 . 5. 2X.F9 .3  > 

0863 

IFC . 5-L£NS>2080. 2080. 2100 

0864 

2080 

IFCS-SLFROI1  )2100.2090.2090 

0865 

2090 

IFC  SLUPT0-S>2100.21 10.21 10 

0866 

2100 

TAe«ARTA3 

0867 

IFCFT.lt. 0.5)  GO  TO  1110 

0868 

IFC AUTON . GT . 0 . 5>  GO  TO  2250 

0869 

GO  TO  1240 

0870 

21 10 

WRITEC6.2120>  T 

0871 

21'20 

F0R«ATC/.9X.*Ta'.2X.F15.6> 

0872 

URITEC6.2130>  XH 

0873 

2130 

F0RHATC9X,*XH**.2X. F15.6  > 

0874 

URITE(6.2140>  XG 

0875 

2140 

F0RHATC9X,*XG-*.2X,F15.6> 

0376 

WRITEC6.2150>  XGH 

0877 

2150 

F0RNATC9X. ’XGN-’. 2X.F15 . 6> 

0879 

UR1TEC6.2160>  YGH 

0879 

2160 

F0RWATC9X.' YGH-* .2X.F15 . 6> 

0880 

URITEC6.2170>  UXF 

0881 

2170 

F0RNATC9X,*UXF»*.2X.F15.6> 

0882 

HRITEC6.2180>  UYF 

0883 

2180 

F0RMATC9X.'UYF»* .2X,F15 .6) 

0  88  4 

HRITEC 6.2190>  UXG 

0985 

2190 

F0RHATC9X. 'UXQ-' ,2X.F15 .6> 

0986 

URI TEC  6. 2200)  UYG 

0887 

2200 

F0RHArC9X. •UYG-' .2X.F15 .6> 

0888 

URITEC6. 2210>  DXG 

0889 

2210 

F0RMATC9X. 'DXG-* .2X.F15 .6> 

0890 

UR1TEC6.2220>  DXGN 

0991 

2220 

F0RNATC9X. ‘OXCN**. 2X. F15 .6 ) 

0892 

UR1TEC6.2230>  DYG 

0893 

2230 

F0RNATC9X.‘DYC»*.2X.F15 .6> 

0894 

UR1TEC6.2240>  REYHN 

0895 

2240 

FORNATCBX. 'REYNH-* , 2X , F 1 5 . 6,  // > 

0896 

GO  TO  2100 

0897 

2250 

BIFCOR-CORTE-ACORTE 

0898 

A8»lC0«A8SCi>IFC0R> 

I  N  T  A 
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0899 

t 

IF(  .  5- ABD  ICO  >2260.- 2260,  1040 

0900 

2260 

IF( D  I FCOR  )227 0,2270, 2280 

090  1 

2270 

YGH03-ftVGHO 

0902 

VGN04-YGH0 

0903 

MRl TE< 6, 2271 >  YGH03.yGH04 

0904 

2271 

FOR«AT</ , 6X, ■ YGH03-*, 2X, F7 . 4 , 2X , ■ YCN04 • ■ , 2X , F7 .4 , /) 

0905 

GO  TO  1040 

0906 

2280 

YGHOl-AYGNO 

0907 

YGH02-YGN0 

0908 

WRI TE< 6, 2281 >  Y6N01,YGK02 

0909 

2281 

FORHATC/ , 6X, ■ YGNOl**, 2X, F7 . 4 , 2X , ■ YGN02- * ,  2X , F7 . 4 , / )  ' 

0910 

GO  TO  1040 

091  1 

2290 

IF<  AUTOff  GT  .  0 . 5)  GO  TO  2310 

0912 

GO  TO  1052 

0913 

2299 

8RI T£< 1, 2300) 

0914 

2300 

F0RflAT<2X, •ATEMTION!  YGHO  OUT  OF  THE’,/, 

0915 

*•  HORIZONTAL  LIMITS  OF  THE  GRAPHIC*,/, 

0916 

*•  00  YOU  MAHT  PRINTED  RESULTS  ?*,/, 

0917 

*•  AHSMER!  YES-1.  ,  H0-0.*,/> 

0919 

READ< 1 ,*  )  PR 

0919 

IF<  PR. GT . 0 . 5  >  GO  TO  1030 

0920 

GO  TO  1240 

0921 

2310 

PP-0. 

0922 

C 

FUTURE  EXPANSION  BY  HERE. 

0923 

2320 

HR!TE< 12 >-l ,0, 9999,5000 

0924 

STOP 

0925 

END 

0926 

BLOCKDATA 

0927 

COHHOH/COH2/X< 100), Y(  100 ) 

0928 

END 

0929 

SUBROUTINE  SU D YG< S , PO 1 , P02 , P03, P04 , P05 , P06 , P07 , 

0930 

•  P08 , P0  9, P  10 , PI  1 , COT, R, S AOt , SA02, SA03,SA04 , SA05 ) 

093  1 

C0HH0H/C0«2.<'X(  100),  Y(  100  ) 

0932 

IF(S)2,1,1 

0933 

1 

IF( . 000001-3)5,5,3 

0934 

2 

IF( . 000001+8)5,5,4 

0935 

3 

S-. 000001 

0936 

GO  TO  5 

0937 

4 

S--  .  00000 1 

0938 

5 

E-S 

0939 

6 

DO  10  1-1,2 

0940 

SSP02-E+P02 

0941 

CALL  SIG<  P03, POl ,SSP02, Pll ,SIG1 ,SIG2,SIG3,SIG4  > 

0942 

VX-<SIGl-SIG3)*P04/3. 141592 

0943 

VY-<SIG2-SIG4  )^P04/3. 14  1592 

0944 

SA04-VX 

0945 

SA05-VY 

0946 

BH-PlO+2 . •<  P04-P06 ) 

0947 

CF-P10*<  P08-2 . *P06  +  VX  > 

0948 

BMP-ABS<  BH**2 . -4 .-CF) 

0949 

8A0l-<BH-BMP**.5>/2. 

0950 

SA0  2  — P05/(  2.  *P04-2  .*P06-SA01  ) 

0951 

SA03-<P09-2 .*P07vVY)*SA0l/<P08-2. *P06+VX> 

0952 

C-SA02*<2.*P07+SA03) 

0953 

ERREL-100 . •ABS<C/E-1  .  > 

0954 

1F<  ERREL-C0T>19,7,7 

0959 

7 

GO  TO  <8, 9), I 

0956 

8 

Cl-C 

0957 

E-C 

0959 

GO  TO  10 

Vim  w 
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* 


I 


1 

0959 

9 

C2-C 

0960 

10 

CONTINUE 

0961 

DCIS-CI-S 

0962 

0C2S=C2-S 

0963 

0C21-C2-C1 

0964 

I>8BS-ABS<  DC21  >-ABS<  DClS  > 

0963 

IF( DC1S>1 1. 11 . 14 

0966 

1 1 

IF( DC21 >12i 12i 13 

0967 

12 

IF( DABS)18. 18. 17 

0968 

13 

IF( 0C2S>ld. 18. 17 

0969 

14 

IF< DC21 >15. IS. 16 

0970 

15 

IF(0C2S>17. 17. 18 

0971 

16 

IF(  l>ABS>18.  18.  17 

0972 

17 

E»S-0C1S**2 . /DC21 

0973 

GO  TO  6 

0974 

18 

E-C2*0C21**2. /DClS 

0973 

60  TO  6 

0976 

19 

R»E 

0977 

RETURN 

0973 

END 

0979 

SUBROUTINE  S1G<NT.EF13.EF14.CMSSF. 

0980 

*SIG1F.SIG2F.SIG3F.SI64F> 

0931 

C0««0H/C0N2/X< 100>.Y( 100) 

0982 

SIG1F»0. 

0983 

SIG2F»0. 

0  93  4 

SIG3F-0. 

0985 

SI64F-0. 

0986 

SC-CMSSF 

0987 

DO  9  J>1 . NT 

0988 

J  1>J 

0989 

J2-J+1 

0990 

XJl -X<  J1 > 

0991 

YJ1»Y<  J1  ) 

0992 

XJ2»X( J2  ) 

0993 

YJ2«Y( J2  ) 

0994 

PVECT»<XJ2-XJ 1  )*<EF14-YJl>-(YJ2-YJl>*<EF13-Xai  ) 

0995 

IFtPVECT  >1. 1.2 

0996 

1 

SC»1. 

0997 

2 

CALL  STER<EF13.EF14.XJ1.YJ1.XJ2,YJ2.TERX.TERY> 

0998 

IF(XJ2-XJ 1 >3. 7.4 

0999 

3 

IF<<YJ2-YJ1 >7<XJ2-XJ1 >>5,8,6 

1000 

4 

IF<  <  YJ2-VJ1  >7<  XJ2-XJ1 >>6.8,5 

1001 

5 

SIGlF-SIGlF-*-TERX*SC 

1  002 

SIG2F«3IG2F*TERY*SC 

1003 

GO  TO  8 

1004 

6 

SIG3F-SIG3F^TERX*SC 

1005 

SIG4F«SIG4F+TERY*SC 

1006 

GO  TO  8 

1007 

7 

CONTINUE 

1  008 

1F( YJ2-YJ 1 >6. 8,5 

1  009 

8 

SC-CNSSF 

1010 

9 

CONTINUE 

101  1 

RETURN 

1012 

END 

1013 

SUBROUTINE  ST ER( EFX . EFY . EFXl . EF Y1 . EFX2 , EF Y2 . SFX. SFY  > 

1014 

XJ-<  EFX1+EFX2  >/2 . 

1015 

YJ-<  EFYl +EFY2  >/2 . 

1016 

PARX-EFX-XJ 

1017 

PARY-EFY-YJ 

1018 

PARX2-PARX««2 . 

tMIMTMU  li*l  1NT4 
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I 


I 


I 


1<)19 

1 

1 

PARY2-PftRY*»2 - 

1020 

0EH0-PARX2+PftRY2 

1021 

VMOO-ABSC  EFY2-EFY1  > 

1022 

SFX»PARXfYHOI>/DEHO 

1023 

SFY-PARY+YMOO/OEHO 

1024 

RETURH 

1  025 

ENO 

1  02^ 

SUBROUTINE  CO R TA< X I , Y 1 , X 2 . Y2 / X3 , Y 3 / X4 / Y4 , XC / YC  ) 

1027 

IF<  XI . Ea  X2  )  GO  TO  1 

1023 

IF( X3. EQ . X4  )  GO  TO  2 

1029 

TERHl-<  Y1-Y2>/'<X1-X2> 

1030 

TERM2-< Y4-Y3)7<X4-X3> 

1031 

TER113«TERM1*K2-TERH2*X3 

1032 

TERfl4-TERHl-TERI12 

1033 

XC=»<  TERH3  +  Y3-Y2)7TERI14 

1034 

YC-TER«1*<XC-X2>+Y2 

103  5 

GO  TO  3 

1036 

1 

XC«Xl 

1037 

YC«< Y4-Y3 )*<XC-X3>/<X4-X3>+Y3 

1038 

GO  TO  3 

1039 

2 

XC>X4 

1040 

YC«<XC-X1 )*<Y2-Y1>/<X2-X1  >  +  Yl 

1041 

3 

CONTINUE 

1042 

RETURH 

1043 

EHO 

1044 

SUBROUTINE  DR AH< REF , XF, YF, NXPF , NYPF > 

1045 

NXPF-500 . •<  4. 5f7 . *XF7REF  ) 

1046 

NYPF«500.*<10.+7  *YF/REF) 

1047 

RETURH 

1048 

END 

1049 

SUBROUTINE  OR AW2( RE F , XDF , Y XDF / XF , YF / NX PF , NYPF > 

1050 

NXPF-500 . *(  4 . 5  +  7 . ♦<  XDF  +  XF )/REF> 

1051 

NYPF-500  ,  10  . +7  .  *<  YXDF  +  YF  )7REF  > 

1052 

RETURH  , 

105  3 

END 

1054 

SUBROUTINE  XYUL< XF , YCF, YTF , S I HHU. 

1055 

*C08NU/XUF,YUF,XLF, YLF> 

1056 

YUF-YCF+YTF*COSHU 

1057 

XUF-XF-YTF+SIHHU 

1058 

YLF«YCF-YTF*COSNU 

1059 

XLF-XF+YTF*SIHHU 

1060 

RETURH 

106! 

1062 

* 

END 

I  N  T  A 
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SUMMARY 

The  ixif ormation  presented  in  this  Pinal  Scien¬ 
tific  Report  can  be  sinmnarized  as  follows: 

1'.-  One  Mathematical  Model  has  been  established, 

determining  the  absolute  motion  of  the  droplets 
composing  a  cloud  when  an  infini'^^ely  long  body, 
whose  arbitrary  cross  section  is  defined  by  up 
to  100  points,  moves  inside  it,  normally  to  its 
infinite  span,  and  that,  in  addition,  determines 
the  trajectories  described  by  these  droplets, 
with  respect  to  the  profile  limiting  the  cross 
section. 

One  parameter,  nominated  NSSP,  has  been 
included  in  the  program  to  take  into  account 
the  possibility  of  attenuate  the  effect,  in  the 
calculation  of  the  air  velocity,  of  the  defining 
segments  that  are  non  sighted  from  the  point 
where  the  air  velocitv  is  being  calculated. 

The  Mathematical  Model  includes  the  ex¬ 
pression  that  determines  the  ice  growing  in  the 
casse  P.P.sl  . 

2.-  One  FORTRAN  Program  ’nas  been  developed  in  base 
to  thah  Mathematical  Model.  This  Program,  nomi¬ 
nated  AIP07  (Aircraft  Icing  Processes,  NO. 7), 
permits  in  both, Manual  auid  Automatic  modes,  to 
obtain  the  trajectories  of  the  droplets  with 
respect  to  the  profile,  presented  in  a  norma¬ 
lized  pattern  at  any  desired  angle  of  attack.- 
The  values  of  13  different  variables  implicated 
in  the  calculations  are  printed,  if  so  desired, 
step  by  step,  along  each  trajectory. 

'When  starting,  the  values  of  up  to  24  para¬ 
meters  that  de-harmine  the  process  and  displav 
of  results  are  to  be  given  to  the  computer  by 
conversation  with  it  using  the  adequate  peri¬ 
pheric. 

As  a  result,  the  AIP07  Program  presents  a 
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great  versatility. 

Iix  Manual  Mode  the  initial  ordenate  of  each 
tra.iectory  is  to  be  selected  by  the  operator. 

On  the  contrary,  in  Automatic  Mode  aJ.1  the 
trajectories  are  executed  automatically  by 
the  computer  according  to  the  data  given  by 
the  operator  before  running  -the  program. 

3.-  Ten  s'ets  of  trajectories,  generated  by  the 
NACA  2409  profile,  are  presented  in  order  to 
show  some  of  the  huge  number  of  possibilities, 
concerning  operation  and  display,  yielded  by 
the  AIP07  FORTRAN  Program. 

The  value  of  the  NSSP  parameter,  ranging 
from  zero  to  one  can  be  selected  for  each  tra>- 
jectory,  or  set  of  trajectories,  and  it  is 
prinxed  at  the  beginning  of  the  tabulation, 
following  the  initial  ordenate  of  each  trajec¬ 
tory. 

The  sheet  that  presents  each  set  of  drawn 
trajectories  is  followed  by  the  printed  sheets 
contSLining  the  corresponding  data  and  results. 
These  results  can  be  the  tabulation  (if  desired) 
step  by  step  of  the  values  of  the  i3  mentioned 
variables  and  information  concerning  if  the 
droplet  hits  the  surface  of  the  profile  or  if 
the  dot  that,  on  the  paper,  represents  the 
droplet  reaches  the  left  hand  side  of  the 
drawing. 
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CONCLUSIONS 


1 The  information  presented  in  this  Pinal 
Scientific  Report  leads  to  the  conclusion  that 
the  Primary  Aim  of  the  proposal  denominated 
"Aircraft  Icing  Processes",  submitted  to  the 
Eiaropean  Office  of  Aerospace  Research  and 
Development  in  February  1933.,  producing 
Graint .NO.AFOSR  83-0340,  has  been  reached. 

2.-  The  AIP07  FORTRAN  Program,  elaborated  under 
that  Proposal,  has  provisions  to  be  expauided 
in  order  to  reach  ne'.7  targets. 

Torre  ,jdn  de  Ardoz,  7  November  1984 
The  Principal  Investigator 
-Srnesto  Montiel  3.- 
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